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It u usually assumed that the iron in iron-contaminated, boron-doped silicon exists as FcB pairs. The iron can be cycled between 
its interstitial (Fe.) and paired (FeB) states with the total density F* + FcB remaining constant. We have discovered that iron 
can also exist in other paired states, which we believe to be Fc-vacancy or Fe-implant damage pairs. When these pairs arc 
destroyed and subsequently when FeB pairs form, we have observed an increased density. We refer to the excess iron f AFeB> 
after re-formation as "latent" iron. • ; 

KEYWORDS: silicon, iron, Impurities, deep-level transient spectroscopy, Fo-B 



Iron is one of the most technologically important impu- 
rities in silicon. Recent reviews detail many of the effects 
of Fe in Si and provide numerous references to previous 
publications. 1 * 25 The International Technology Roadman for 
Semiconductors requires iron densities below 10 in cm~* for 
future Si processing. ^ Detection of such low iron rlftn.ciH^ 
within the Si wafers is possible by several characterization 
techniques, e.g. 9 surface photovohage (SPV), photoconduc- 
tance decay (PCD), and deep-level transient spectroscopy 
(DLTS). 4> The former two rely on minority carrier diffusion 
length and recombination lifetime measurements using the ki- 
netics of Fe-B dissociation in boron-doped silicon as a finger- 
print,^ while the latter is a direct measurement. Since Fe in 
boron-doped silicon exists as either Fe-B pairs or as intersti- 
tial iron, Fej, DLTS allows a direct measurement of iron in 
either one of these states. In this paper we show that these 
characterization techniques can underestimate the Fe density 
due to a third iron state— the latent iron. The error can be as 
high as 100%, which can be critical for today's silicon, since 
such device parameters as diode leakage current, gate oxide 
integrity, and carrier lifetime depend strongly on iron density. 

We have discovered a new and previously unreported effect 
associated with interstitial iron. This is an increase of the iron 
density in the wafer after annealing, as illustrated in Fig. 1. 
This figure is a set of DLTS spectra of a Fe-contaminated Si 
wafer. Trace (a) is the original spectrum aAer wafer room 
temperature storage for two years. It clearly shows the well- 
known FeB peak at T = 55.3 K and essentially no interstitial 
iron at T = 253 K. After a FeB dissociation anneal at 1 80°C 
for 30 s, the FeB peak disappears and the Fe; peak appears 
at T st 253 K, as shown by trace (b). The magnitude of the 
two peaks is approximately the same, as it should be since 
Fe-B pairs are changed into Fe ; and substitutional boron B e . 
Holding the sample at 80°C for ) 5 h followed by room tem- 
perature storage for 5 days, causes Fcj and B, to change back 
into FeB pairs, shown by trace (c). The change from Fcj into 
FeB pairs is well known. What has not been reported before 
is the significant increase of the FcB density from the initial 
2 x 10 ,2 cm- 3 (trace (a)) to the final 4.5 x IO ,2 cm" 3 (trace 
(c)). We call this excess iron "latent iron'\ because it was not 
ongtnally detected as a FeB pair density. What is the source 
of this excess iron? 

We used double side-polished float-zone (FZ), p-type sir- 
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Fig. 1. DLTS spectra of(a)ihc original sample, (b) after 1 80*C, 30 s disso- 
ciation anneal, and (c) after StTO 1 5 h and room temperature 5 day storage. 



icon wafers doped with boron to 2-1 x I0 u cm~ 3 . Iron 
was introduced by thermal annealing at 900°C for 2 hours 
following iron implantation into the back surface with a 
dose of I0 n cm -2 and energy of lOOkeV. DLTS mea- 
surements on the front wafer surface yielded an iron den- 
sity of 2 x 10 l2 cm- 3 , Schottky diodes were fabri- 
cated by aluminum-evaporated contacts on the front sur- 
face and a large-area silver paste contact on the back sur- 
face. The samples were characterized with; current-voltage, 
capacitance-voltage, and DLTS measurements. 

Capacitance-voltage data showed little change before and 
after dissociation anneal. This is expected since the effective 
doping density change, due to boron becoming inactive in the 
FeB state and active in the unpaired FeB state, is very small 
because the iron density (- I0 12 cm~ 3 ) is much lower than the 
boron density 10 U cm" 3 ). Schottky diode current-voltage 
measurements showed a substantial current and diode ideality 
factor change upon dissociation anneal, possibly associated 
with changes at the metal/Si interface. 

The key DLTS results are shown in Figs. 2(a) and 2(b). 
The samples were initially measured after storage at room 
temperature for 2 years. This is the "As-is" curve in both 
figures. Then they were annealed at IS0°C/30s to dissoci- 
ate the Fe-B pairs. The FeB peak vanished and the Fe* peak 
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(b) 

Fig. 2. DLTS spectra for "As-is", after 1 80°C/30 s dissociation anneal and 

(a) 80° C storage for various times followed by room temperarure/5 days, 

(b) roam temperature storage for various times. 



appeared. The wafer in Fig. 2(a) was then annealed at 80°C 
for various times up to 1 5 h and DLTS spectra were recorded. 
Then the sample was held at room temperature for five days 
and a last DLTS spectrum was taken. The Fe* peak has essen- 
tially disappeared and the FeB peak has reappeared, but it is 
significantly higheT than the original. The sample in Fig. 2(b) 
underwent a dissociation anneal as that in 2(a), but was held 
at room temperature thereafter, instead of at 80°C, and DLTS 
spectra were recorded. The results are qualitatively similar 
to those in Fig. 2(a). Note that the magnitude of the "As-is" 
Fe-B peak and the Fcj peak immediately after the dissociation 
anneal are about the same. It is only after an extended time at 
room temperature or 80°C that the latent iron appears in the 
new Fe-B peak. The initial 180°C/30s dissociation anneal is 
insufficient to generate the latent iron. 

In both cases we note the FeB density to be higher after this 
anneal sequence than at the start. This "latent" iron, denoted 
as AAV,) B in Fig. 1, is higher for the 80°C anneal than for the 
room temperature anneal case. Wc express this sequence as 

(FeB) 0 -> (Fei) + + (B,)~ a(FcB)° 



B.-D. Choj tff o/, 

with a between 1 and 3. For the T = 80°C anneal, a 2. 
The temporal behavior is shown in Fig. 3 for the 80°C anneal 
case. The FeB density starts out at 2 x 10 ,2 cm" J . drops to 
approximately zero after the 180°C/30s anneal, then gradu- 
ally increases and begins to saturate at 5.4 x 10 4 s (15 h), the 
limit of our measurement. After storing the sample at room 
temperature for 5 days, 7/ F cB increases further. A# FaB . de- 
fined as the final density minus the initial density, is plotted 
versus anneal temperature in Fig. 4. It remains reasonably 
constant, then increases rapidly at around 70°C and appears 
to saturate around 80~90 V C. For each datum point in Fig. 4, 
AJV> cB is defined as in Fig. 3. At this point we do not un- 
derstand the rapid A^ eB increase from 70°C to 80°C. From 
Fig. 3 it appears the additional 5 day room temperature anneal 
contributes significantly to this increase. 

It is generally assumed that all of the iron in boron-doped 
Si is completed with boron after storage at room temperature 
for a few days or longer/* - ** However, in addition to boron, 
there may be other possible sinks for iron. Two of these 
are growth-induced thermal vacancies, and ion implantation- 
induced damage. Since vacancies can be negatively charged 
in silicon, they may have Coulombic interaction with posi- 
tively charged Fcj. Furthermore, Fe|**-V pairing may be fa- 
cilitated by the difference in covalent tetrahedral radii of Fe 
and Si atoms. A simple calculation shows that the density of 
thermal vacancies after a 900*0 anneal is sufficiently high to 
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Fig. 3. Nt cii as a function of time following various anneals. 
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Fig. 4. AJVfcB as a function of anneal tempera ture. The time for each tem- 
perature was 1 5 h ai each temperature + 5doys ot room temperature. 
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bond with 2 x 10 12 cm" 3 iron atoms. 

Since Fe was introduced in our samples by ion implanta- 
tion at 100 kcV, it is possible that implantation-induced dam- 
age occurred in the wafers. This damage could consist , of 
silicon vacancies and interstitial*. The majority of these point 
defects may be annihilated during the 900°C anneal, but some 
of them may cluster in the form of microdefects or dislocation 
loops to serve as sinks for iron. The implant created vacancies 
are generated at the side opposite to the DLTS measurements. 
The vacancy diffusion coefficient is difficult to determine, but 
has been given as 1CT 5 cm 2 /s at 1 100°C with an activation 
energy of 0.35 to 0.45eV 9 l0) These values lead to a diffu- 
sion coefficient of about 5 x 10~ 6 cm 2 /s al 900°C, giving a 
(2Dt) ,/2 diffusion length of 2700 *im, more than sufficient to 
diffuse through our 700 fim thick wafers. The iron was intro- 
duced by thermal annealing at 900°C for 2 hours following 
iron implantation. 

We envisage that during the 1 80 P C anneal, some of the iron 
atoms may be released from the two types of sinks and be- 
come available to form FeB pairs. Even though we are not 
able to distinguish between the two possible sinks, we feel 
that ion implantation-induced damage may dominate as a trap 
for iron atoms because it is more extensive. Interestingly, 
these pairs do not appear on the DLTS plots over the tem- 
perature range of 40-300 K, because they cannot be detected 
by DLTS or their energy level is such that it is not accessible 
in the temperature range scanned in the DLTS measurements. 
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To summarize, we have observed iron in boron-doped Si 
that we believe to be paired with defects rather than with 
boron. First, we determine the density of FeB pairs, then 
dissociate these pairs by a 180°C/30s anneal and store the 
samples sufficiently long to reform FeB pairs. The reformed 
FeB pair density is higher by a factor of two or more than 
the initial density. We believe the Fe-dcfcct pairs to be either 
Fe-vacancy or Fe-implant damage pairs. 

The research leading to this paper was partially funded 
by the Silicon Wafer Engineering and Defect Science 
Consortium (SiWEDS) (Intel, Komatsu Electronic Metals, 
LG Siltron, MEMC Electronic Materials, Mitsubishi Silicon] 
Okmetic, Sumitomo Siiix Silicon, Texas Instruments, and 
Wacker Siltronic Corp.). 
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Abstract This article is the first in a series of two reviews 
on the properties of iron in silicon. It offers a comprehensive 
summary of the current state of understanding of fundamen- 
tal physical properties of iron and its complexes in silicon. 
The first section of this review discusses the position of iron 
in the silicon lattice and the electrical properties of interstitial 
iron. Updated expressions for the solubility and the diffusiv- 
ity of iron in silicon are presented, and possible explanations 
for conflicting experimental data obtained by different groups 
are discussed. The second section of the article considers the 
electrical and the structural properties of complexes of inter- 
stitial iron with shallow acceptors (boron, aluminum, indium, 
gallium, and thallium), shallow donors (phosphorus and ar- 
senic) and other impurities (gold, silver, platinum, palladium, 
zinc, sulfur, oxygen, carbon, and hydrogen). Special attention 
is paid to the kinetics of iron pairing with shallow accep- 
tors, the dissociation of these pairs, and the metastability of 
iron-acceptor pairs. The parameters of iron-related defects in 
silicon are summarized in tables that include more than 30 
complexes of iron as detected by electron paramagnetic res- 
onance (EPR) and almost 20 energy levels in the band gap 
associated with iron. The data presented in this review il- 
lustrate the enormous cornplexing activity of iron, which is 
attributed to the partial or complete (depending on the tem- 
perature and the conductivity type) ionization of iron as well 
as the high diffusivity of iron in silicon. It is shown that 
studies of iron in silicon require exceptional cleanliness of 
experimental facilities and highly reproducible diffusion and 
temperature ramping (quenching) procedures. Properties of 
iron that are not yet completely understood and need further 
research are outlined. 

PACS: 71.55;Cn; 76.30.Fc; 76.80.+y; 66.75.+g; 78.30.Am; 
66.30Jt; 78.55.Ap; 81.05.Cy; 61-.72.-y 



Iron is one of the most ubiquitous and detrimental metal im- 
purities in silicon. For example, in metal-on-insulator (MOS) 
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structures, iron can precipitate at the Si-SiC>2 interface, thus 
locally thinning the oxide and increasing surface roughness. 
This causes premature electric field breakdown of oxides and 
degrades the gate oxide integrity (GOI) [1-7]. In the deple- 
tion regions of MOS devices, iron acts as an effective minor- 
ity carrier generation site [7,8], increasing the dark current of 
charge-coupled devices (CCD) [9] and causing refresh fail- 
ures in dynamic random access memory (DRAM) chips [6]. 
The generationnrecombination centers associated with dis- 
solved iron and its complexes increase leakage currents in any 
reverse-biased junction, thus increasing power consumption 
and heat production [10-15]. In crystalline and polycrys- 
talline photovoltaic devices, iron contamination results in the 
creation of recombination centers which reduce the minority 
carrier diffusion length and consequently the solar cell ef- 
ficiency [16]. As device dimensions continue to shrink, the 
device yield is becoming ever more sensitive to defects and 
impurities [17,18]. Since the ultrapure technology that en- 
ables one to reduce iron surface contamination to the level of 
below 10 10 to 10 n cm" 2 is extremely expensive, it has be- 
come increasingly important to understand the mechanisms 
of iron contamination, the detrimental role of iron in silicon 
devices, and the tolerable limits of iron contamination for 
each particular technological process. This demand has stim- 
ulated great interest in the physics of iron in silicon. 

The major obstacle in establishing a coherent experimen- 
tal picture of 3d transition metals in silicon, and of iron in 
particular, arises from the high diffusivity of these elements. 
Iron remains mobile at room temperature and can diffuse 
quickly at temperatures above 100 °C. Since interstitial iron 
is positively charged in p-Si at room and slighdy elevated 
temperatures, it readily forms pairs with negatively charged 
defects, such as shallow acceptors. The high reactivity of iron 
was confirmed by electron paramagnetic resonance (EPR) 
studies, which detected more than 30 complexes of iron, and 
by electrical measurements, which revealed about 20 deep 
levels associated with iron (see the discussion below and ta- 
bles in Appendices A, B). 

Iron-related defects in silicon have been intensively 
studied for more than thirty years. Substantial progress has 
been made in the last ten to fifteen years. Unfortunately, no 
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comprehensive review of the properties of iron in silicon has 
been published since the 1983 review of Weber [19] (and 
some authors even reference the early reviews of Milnes and 
Chen [20,21]). The recently published book by Graff [22] 
is the only exception, but since it covers the whole scope of 
metal impurities in silicon device fabrication, it is limited to 
a rather short discussion for each particular metal and does 
not offer an in-depth insight into the physics of iron in silicon. 

Our review is intended to cover the gap between the 1983 
review of Weber [19] and the current state of knowledge by 
giving a comprehensive up-to-date analysis of electrical prop- 
erties of iron in silicon, its structure, complexing behavior, 
and recombination properties. We intentionally left out of the 
scope of this article such important issues as the degrada- 
tion of devices due to iron contamination, the sources of iron 
contamination in modern silicon technology, methods of de- 
tection of trace concentrations of iron, and the gettering of 
iron. These will be covered in a separate review, which will be 
published in one of the forthcoming issues of this journal. 



1 Isolated iron as a defect in silicon 

LI Lattice site of iron in silicon 

The donor behavior of iron (see Sect 1 .2,2), its high diffusion 
coefficient, and the agreement of the electron paramagnetic 
resonance (EPR) of iron with the model for EPR resonances 
of interstitial 3d metals enabled Woodbury and Ludwig [23] 
to conclude that Fe diffuses interstitially. The interstitial lat- 
tice site of iron was confirmed by Greulich-Weber et al. [24] 
based on the agreement between the experimentally deter- 
mined parameters of superhyperfine interactions of iron with 
silicon ligands measured by electron nuclear double reson- 
ance (ENDOR) measurements, and theoretically calculated 
parameters for interstitial iron (Fej). The isomer shift and 
the broadening of the Mossbauer lines observed by Schwal- 
bach [25] are also in line with the interstitial lattice site of iron 
in silicon. 

There is much speculation in the literature as to what frac- 
tion of the total dissolved iron concentration is substitutional 
iron (Fes). For instance, Collins and Carlson [26] observed 
a "slow" and "fast" diffusing 59 Fe species, and assigned 
them to an interstitially and substitutionally (via vacancy or 
interstitial-assisted mechanism) diffusing iron, respectively. 
However, later Stojic et al, [27] demonstrated that the "slow- 
diffusing" component actually arises from iron that precip- 
itated at the sample surface. To the best of our knowledge, 
there is no direct experimental evidence for the existence of 
substitutional iron in a non-distorted configuration. More than 
30 iron-related complexes have been detected by EPR, but no 
iron-related resonance has been found with a g-tensor com- 
patible with the model of Ludwig and Woodbury [28] for 
substitutional metals. The only complex that can be inter- 
preted as substitutional iron in a distorted configuration was 
detected by Muller et al. [29] in irradiated iron-doped sili- 
con samples. Hyperfine splitting of the spectrum due to 57 Fe 
clearly showed the presence of one iron nucleus in the center. 
Since the second component of the complex in irradiated sil- 
icon is most probably a vacancy, it was suggested that this is 
a singly positively charged interstitial iron-vacancy complex, 
(FeiV) + . Because the center has (111) axial symmetry, the 



iron atom can not be exactly at the substitutional position but 
has to be shifted along a trigonal axis towards the tetrahedral 
interstitial position [29] (see Appendix A). 

He et al. [30] reported that the co-diffusion of Fe and 
Cu in p-type silicon substantially decreases the concentra- 
tion of interstitial iron detected after diffusion, although no 
new energy levels associated with FeCu pairs were found in 
the lower half of the band gap. It was suggested that the de- 
crease of the concentration of Fej is due to the transformation 
of interstitial iron into substitutional iron, caused by the injec- 
tion of vacancies by diffusing and precipitating copper [31]. 
While this model agrees with the observations of Woodbury 
and Ludwig [32] who reported that co-diffusion of Cu with 
manganese or chromium (or co-diffusion of silver with Mn) 
resulted in the transfer of interstitial Mn and Cr onto substi- 
tutional sites, it is not clear why diffusion of copper should 
generate vacancies since it is known that copper silicide cre- 
ates a strong compressive lattice strain and is thus expected to 
generate interstitials [33,34]. 

Gilles et al. [35] found a substantial increase (by three 
orders of magnitude) of iron solubility in n + -Si with a phos- 
phorus concentration 1 x 10 20 cm~ 3 along with a decrease in 
the iron diffiisivity by almost four orders of magnitude. This 
was interpreted as an increase of either the fraction of sub- 
stitutional iron due to heavy n-type doping, or the fraction 
of iron paired with phosphorus [35]. A substantial increase 
of substitutional solubility in n + -Si can indeed be expected 
if substitutional iron forms acceptor level(s) in the band gap. 
Hall and Racette [36], following Shockley and Moll [37], 
showed that the substitutional solubility, S s , of transition 
metals in extrinsic silicon varies in the nondegenerate case as 
S s = Sso x [n/n\Yy where is the substitutional solubility 
in intrinsic silicon, n and n\ are the free electron concentra- 
tions in doped and undoped samples, respectively, and r is the 
multiplicity of the acceptor level of the metal. However, there 
is no experimental evidence in favor of existence of an accep- 
tor level of substitutional iron. In fact, theoretical calculations 
using spin-unrestricted density-functional theory with the lin- 
ear muffin-tin-orbital Green's function method [38,39] as 
well as self-consistent local-density-functional calculations 
using ab initio nonlocal pseudopotentials [40,41] predicted 
that although neutral substitutional iron should be stable in 
the silicon lattice, it should not have any electrical levels lev- 
els in the band gap. A similar conclusion was made from 
the analysis of Mossbauer spectroscopy data (see, for ex- 
ample, [42-44] and references therein), which revealed one 
line in the spectra that was attributed to interstitial iron and 
another one attributed to substitutional iron. Weyer et al. [42] 
pointed out that the same isomer shift determined for the 
line attributed to substitutional iron in p 9 «, and n + silicon 
is consistent with the above mentioned ([38—41]) theoretical 
conclusion that Fes has no band gap state [38-41]. Note that 
Mossbauer studies have to use ion implantation to achieve 
a sufficiently high local concentration of iron, and a substi- 
tutional relocation due to a replacement collision with the 
lattice atoms was suggested as the source of substitutional 
fraction of iron [42]. 

Thus, there are no convincing experimental data that sub- 
stitutional iron exists in detectable concentrations in moder- 
ately doped silicon, at least if no vacancies were generated 
by irradiation or implantation. There are also almost no data 
on the properties of iron in n + -Si where the existence of high 
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concentrations of Fes is viable if substitutional iron forms an 
acceptor level For this reason, we will confine ourselves in 
the following discussion only to the properties of interstitial 
iron and its complexes. 

12 Electrical level and microscopic structure of interstitial 
iron in silicon 

12.1 Microscopic structure of interstitial iron: EPR studies. 
A model for the electron structure of the 3d group impuri- 
ties in silicon was developed by Ludwig and Woodbury [28, 
45,46]. According to their model, illustrated in Fig. 1 for the 
case of iron, all 4^-shell electrons are transferred to the 3d 
shell. Then, all N valence electrons (filled circles in Fig. 1) 
of the impurity (atomic configuration 3d m 4s n with N = 
rn + n) fill the levels (indicated by lines in Fig. 1) in such 
a way that a maximum electron spin S is achieved (Hund's 
rule). The tetrahedral symmetry of the crystal field splits the 
tenfold degenerate atomic 3d orbital into a sixfold degenerate 
h orbital and a fourfold degenerate e orbital separated by the 
crystal-field energy (see Fig. 1). Finally, the model of Ludwig 
and Woodbury assumes a sizeable reduction in the spin-orbit 
interaction constant. Such "quenching" of the spin-orbit in- 
teraction is characteristic for all 3d metals. It is explained 
by either hybridization of the transition-metal orbitals with 
the surrounding ligands [47], or by the dynamic Jahn-Teller 
effect [48-50]. 

The EPR spectrum of neutral interstitial iron, Fef , was 
discovered by Ludwig and Woodbury [23,51] and Feher [52]. 
They reported that the spectrum of neutral interstitial iron of 
natural isotopic abundance consists of one broad line with 
a poorly resolved structure, spin S = 1, cubic symmetry, and 
g = 2.07 (see Appendix A). For higher microwave energy, 
a narrow peak appears in the center of the broad EPR spec- 
trum. This spectrum was explained by attributing the broad 
line to the usual AM = ±1 transition and the narrow line to 
the double quantum transition AM = 2 (Af being the pro- 
jection of a total angular momentum of the Fej along the 
direction of the magnetic field). For magnetic field jdirec- 
tion parallel to (100) and high microwave power levels, the 
central peak measured in the dispersion mode was found to 
show a pronounced structure of at least seven equidistant lines 
which were assigned to the hyperfine interaction with 29 Si nu- 
clei in two shells of the nearest and next-nearest neighbors 
of the interstitial iron atom [53]. The broad line could be re- 
solved into two components by applying a uniaxial stress to 
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Flg,l. The electronic structure of ionized and neutral interstitial iron in 
silicon according to the model of Ludwig and Woodbury (after [28,46]) 



the sample [23,54]. This can be easily understood since the 
zero-field splitting of this 5=1 center is negligible in full 
tetrahedral symmetry and the M = 1 0 and M = — 1 <+ 0 
transitions coincide. This degeneracy is lifted by applying 
a uniaxial stress [54,55]. 

The isotropic EPR line of Fej** was first observed in p-Si 
by Woodbury and Ludwig [23] and was identified in their 
next paper [56] (see Appendix A). According to the model 
of Ludwig and Woodbury [28] (see Fig. 1), the ionized in- 
terstitial iron should have three unpaired electrons and a spin 
of 5 = 3/2. However, the effective value of S = 1/2 is also 
frequently used in the literature (for example, [46,57]). Van 
Kooten et al. [58] showed that Fej^ interacts strongly with 
the silicon lattice, as the electron spin distribution can be 
measured over a volume of 98 atoms. From the hyperfine in- 
teractions, they derived that the lower limit of the amount of 
transferred spin to these ligands is 26%. For neutral iron, it 
was found that from 5% to 25% of the spin density is delocal- 
ized on the first six shells of silicon neighbors [24,59]. 



122 Electrical properties of interstitial iron. Although nu- 
merous experimental data, obtained on intentionally and un- 
intentionally iron-contaminated samples [60-74] suggested 
that the level at E r ^ Ey + (0.40± 0.05) eV (for briefness, in 
the following we will use the notation Ey + 0.40 eV instead of 
Et £ v + 0.40 eV) is associated with iron, it took about 20 
years of research starting from the pioneering work of Collins 
and Carlson [26] to prove that this level is indeed the level of 
interstitial iron. The ultimate evidence for this was obtained 
by Feichtinger et al. [75], who diffused iron in n- and p-type 
silicon with different resistivities and plotted the intensity of 
the EPR signal as a function of the calculated Fermi level pos- 
ition (Fig. 2). Clearly, the intensity of the EPR spectrum of 
Fef depends on the fraction of iron that is neutral at a given 
temperature, and thus on the Fermi level position in the sam- 
ple. As follows from Fig. 2, the charge state of Fej changes 
from Fej* - to Fef as the Fermi level crosses the interval from 
Ey + 0.36 eV to Ey + 0.45 eV. A more exact data evaluation, 
which took into account the temperature dependence of the 
intensity of the EPR lines, enabled Feichtinger et al. [75] to 
conclude that the position of the interstitial iron level is Ey + 
(0.375 ±0.01 5) eV. Additional evidence supporting this con- 
clusion include the identical annealing behavior of the Fej 1 
center detected by EPR and that of the trap at about Ey +■ 
0.4 eV [76], the observations that the concentration of this 
level is close to the total iron concentration measured by neu- 
tron activation analysis and EPR [77,78], and that only this 
level was detected after rapid quenching of samples contam- 
inated with iron [79]. 

The donor nature of interstitial iron has already been 
established by the first Hall effect and resistivity measure- 
ments [26,80-82], which indicated a reduction in the density 
of free holes in quenched samples contaminated with iron. 
Collins and Calson [26] and Wunstel et al. [80] concluded 
from their Hall effect measurements that the contamination 
of n-Si with iron does not result in any substantial change in 
the concentration of free electrons, thus indicating that there 
is no acceptor level of iron, at least at concentrations above 
10 14 cm" 3 . Feichtinger [75] provided experimental evidence 
that, within the limits from E c - 0.045 eV. to Ey + 0.045 eV, 
interstitial iron exists in only two charge states: Fef and Fej*'. 
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Fig Z The intensity of the EPR signal associated with Rf versus Fermi 
level position at 95 K (after Feichtinger et al. [75J). The data points 
^ Uf 5"P boron aluminum, or phosphorus-doped silicon with differ- 
ent levels of doping wh,ch enabled the authors of [75] to identify the energy 
level of the interstitial iron ™<«*y 



The temperature dependence of the position of iron level 
in silicon has received increasing attention in the last several 
years due to its importance for modelling of segregation-type 
gettenng in p/p+ epitaxial wafers. The available experimen- 
tal data indicate that the energy levels of iron, as well as 
of manganese and cobalt are strongly temperature depen- 
dent at high temperatures. These data were inferred from 
the analysis of the solubility of these elements as a func- 
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Fig.3. The temperature dependence of the position of the donor level of iron 
in the band gap with respect to the valence band edge. The symbols 'repre- 
sent the experimental data for the iron level position obtained in numerous 
DLTS and Hall effect studies at low temperatures (filled square) and as re- 
ported by McHugo et al. [83.84] {filled circles) and Gilles et al. [35] (open 
triangles). A data point measured by McHugo et al. [83,84] at 1373 K cor- 
responded to the iron level merging with the valence band and is shown 
in this figure by an arrow. The lines represent attempts to fit the tempera- 
ture dependence of the iron level by three different functions: a parabolic 
fit (curve /, dotted line), a step-like function (curve 2, solid line) and a pro- 
portional fit EMT)/E%(T) s const (curve i, dashed line). The temperature 
dependence of the band gap width E g (T) is shown by the thick solid line 
"Ac" This dependence does not take into account an additional band-gap 
shrinkage due to heavy doping 



tion of doping level [35,83,84]. McHugo et al. [83,84] and 
Gilles et al. [35] concluded from their neutron activation an- 
alysis (NAA) and Mossbauer spectroscopy studies, respec- 
tively, that the donor level of iron moves toward the valence 
band as the temperature increases. Their experimental data 
are presented in Fig. 3. Both McHugo et al. [83,84] and Gilles 
et al. [35] concluded that the position of the iron level at 
T = 800 °C coincides with that measured by DLTS and Hall 
effect, whereas at T> 900 °C the iron level dives towards the 
valence band and nearly merges with it at T> 1 100 °C This 
trend follows from both sets of data. McHugo et al. [83,84] 
pointed out that the iron level shifts so closely to the valence 
band at 1100°C that the Boltzmann approximation was no 
longer valid, so he could not calculate the position of the 
iron level at 1 100 °C. The absence of experimental data be- 
tween the temperatures where the DLTS and Hall effect data 
were obtained, and 1000 K, which was the lowest tempera- 
ture of the solubility studies, makes the interpolation of the 
iron level position over the whole temperature range difficult. 
Since in the next sections we will repeatedly make estimates 
of the ionized fraction of iron at different temperatures and 
for various types of doping and doping level, it is import- 
ant to make certain assumptions about the possible shape 
of the temperature dependence of the iron level despite the 
lack of experimental data. One can imagine three substan- 
tially different simple curves for a temperature dependence of 
the iron level in the intermediate temperature range, as indi- 
cated in Fig. 3: a parabolic fit to the experimental data (curve 
1), a step-like function (curve 2) and a 1 smooth curve given 
by (£Fe(T) - E v )/E g (T) = const (curve 3), where E g (T) is 
the temperature-dependent band gap width. In the first case, 
the depth of the iron level with respect to the valence band 
edge initially increases and then decreases at high tempera- 
tures. In the second case the iron level remains constant up to 
T = 1000 K and then dips rapidly toward the valence band. In 
the last case it drifts smoothly downwards in the whole tem- 
perature range. The proportional dependence (curve 3) agrees 
poorly with the experimentally observed trend in iron level 
position at high temperatures. A step-like function (curve 2) 
can describe the experimental data quite well. However, an 
abrupt change of slope in curve 2 would imply some sort 
of phase transition [85] in the FerSi system. There is clearly 
no evidence for such a transition. Thus, the parabolic depen- 
dence is the only curve of the three, which is smooth and fits 
well the experimental data. Although it would be premature 
and speculative to conclude that the temperature dependence 
of the iron level position in the band gap is described by 
this parabolic function, we want to attract attention of the 
scientific community to this possibility. Experimental studies 
at T<z 1000 K are required to determine which of the three 
models provides a better description for the temperature de- 
pendence of iron in silicon. 

The emission of holes from the level of Fej was found 
to be strictly exponential [86]. Despite this, the energy level 
of Fej as determined from deep level transient spectroscopy 
(DLTS) measurements [78,80,87-90] varied in different 
studies from £ v +0.39eV to £ v + 0.45eV. This scatter of 
emission activation energies is caused by the thermally ac- 
tivated hole capture by Fef with an activation energy from 
0.04 eV to 0.05 eV, as reported by Brotherton et al. [91] (<r p = 
5.6 x 10- l6 exp(-0.048eV//:B^), cm 2 ), Wunstel et al. [80] 
(a p = 1.6 x 10~ 16 x exp(-0.043 eV/* B T), cm 2 ), and In- 
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Fig.4. Temperature dependence of the hole capture cross-section of intersti- 
tial iron in silicon according to [66,80,91,92]. The solid line represents a fit 
to all data points, as discussed in the text 

dusekhar et al. [66] (a p = 6.8 x 10" 16 x exp(-0.04 eV/k B T), 
cm 2 ). The data on toe hole capture cross-section of Fej, 
including the data from [90,92,93], are summarized in 
Fig. 4. This figure does not include the data of Gerson 
et al. [60], who reported a hole capture cross-section of 
CT P = 4 x 10~ 18 cm 2 but did not explicitly specify the tempera- 
ture at which it was measured. Graff et al. [78] pointed out 
that the hole capture cross-section of a p = 3 x 10™ 16 cm 2 re- 
ported earlier by Schibli and Milnes [94] agrees with then- 
data, again without specifying the temperature. 

The data points in Fig. 4 are strongly scattered. How- 
ever, the values for the capture barriers reported in [66,80,91] 
are similar, which suggests that the scatter is chiefly in the 
pre-exponential factor. Therefore, we averaged the capture 
barriers reported in [66,80,91] and then determined the pre- 
exponential factor by making a least-squares fit to all data 
points in Fig. 4. This gave us the following expression for the 
hole capture cross-section of interstitial iron: 

Op(Fei) = (3.9±0.5)xKr 16 

/ 0.045 ± 0.005 eV\ . 9 

xexp ( Sr J (cm) * (1) 

It is worth noting that the barrier for hole capture of 0.045 eV 
explains the discrepancy between the energy level of iron at 
£ v + 0.38, obtained from EPR [75] and Hall effect [21,80], 
and the average value of Ey + 0.43 eV, obtained from tem- 
perature dependence of majority carrier emission rate from 
the iron trap measured with DLTS [80,90,91]. This is because 
the evaluation of DLTS data by means of an Arrhenius plot 
does not enable one to distinguish the emission enthalpy and 
the barrier for capture, and instead yields the sum of both 
energies: 0.38eV+0.045eV = 0.425 eV. To the best of our 
knowledge, the nature of this capture barrier was not dis- 
cussed in the literature. 

The optical ionization energy of the level of interstitial 
iron, as determined by optical methods such as photoconduc- 
tivity, photocapacitance, and deep level optical spectroscopy, 



is close to the thermal ionization energy. The reported pho- 
toionization energies for the transitions between the impurity 
level and valence band are in good agreement with each other: 
0.43 eV [95], 0.42 eV [66], 0.39 eV [87], and 0.38 eV [81] 
are close to the thermal ionization energy. The data for tran- 
sitions from the defect level to the conduction band also 
show a relatively small scatter: AE opt = 0.75eV [66,95- 
102];0.76eV [78,87,103]; 0.78 eV [81,104,105]. 

The electron capture cross-section of Fe^ is not studied as 
accurately as the hole capture cross-section. Collins and Carl- 
son [26] studied the influence of iron on the minority carrier 
lifetime in n-Si and /?-Si and found that a n > 1 .5 x 1 0" 15 cm 2 
at room temperature. Lagowski et al. [106] and Zoth et 
al. [107] reported an electron capture cross-section of Fej of 
cr n = 4x 10~ I4 cm 2 at 300 K. 

A number of theoretical calculations for the structure of 
electron orbitals and deep levels associated with interstitial 
iron in silicon have been reported. One of the first calculations 
was done by Roitsin and Firshtein [108], who used a LCAO 
molecular orbital method and predicted the position of the 
t level above the e level, in contradiction to the model of 
Ludwig and Woodbury (see Fig. 1). However, later calcula- 
tions using the spin-restricted scattered-wave Xct method re- 
ported by DeLeo et al. [109] and Watkins et al. [110], and the 
self-consistent local-density Green's-function calculations of 
Zunger et al. [40, 111] confirmed that the 5t 2 level of iron lies 
below the 2e level. It was shown that the 5t 2 and 2e levels rep- 
resent the antibonding states resulting from the interaction of 
the transition metal 3d state with silicon levels near the top of 
the valence band. 

The positions of the Fej level in the. bandgap predicted 
by the theoretical studies mentioned above show a fair agree- 
ment with the experimental data. DeLeo et al. [109,112] 
predicted the Fep + position at E Y + 0.68eV; Zunger et 
al. [40] calculated the position of the Fe? /+ lever at about 
Ey -f 0.53 eV; Beeler et al. [38,39] predicted the Fe? /+ level 
at fs Y +0.25eV. The agreement between experimental data 
and theoretical predictions became better as the theoretical 
models developed and the available computational power of 
the computers increased. Katayama-Yoshida and Zunger [47, 
113-115] used the self-consistent all-electron spin-polarized 
Green *s-fiinction method in the framework of the local-spin- 
density (LSD) formalism and predicted the iron level at Ey + 
0.32 eV, whereas Weihrich and Overhof [116] predicted the 
iron level position at Ey + 0.29 eV from ab initio total energy 
calculations. However, reliable predictions of the energy lev- 
els with an accuracy better than 0.1 eV are still not available 
with today's theoretical approaches. 



1.3 Diffusivity of interstitial iron in silicon 

1.3.1 Charge state of diffusing iron. The diffusion coeffi- 
cients of transition metals in silicon depend on the charge 
state of the metal since the ionic radii of the neutral and 
ionized species are different, and the degree of the elec- 
tronic density overlap of the diffusing atoms with surround- 
ing silicon matrix is different. Migration enthalpies and pre- 
exponential factors of the largest of 3d ions can be astonish- 
ingly accurately calculated as a difference in the elastic ener- 
gies at the tetrahedral and the hexagonal interstitial sites using 
a hard-sphere model, as suggested by Utzig [117]. However, 
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this model fails for Cu, Ni, Co, and Fe, which have the small- 
est ionic radii from the 3d row. For instance, for neutral 
iron, Utzig [117] predicted a diffusion barrier of 0.46 eV and 
a pre-factor Do = 1.34 x 10~ 3 cm 2 s -1 . No prediction could 
be made for Fef since the ionic radius of positively charged 
iron is too small to distort the lattice. Obviously, the diffusion 
barrier of interstitial iron is determined primarily by a strong 
interaction of iron valence electrons with up to a hundred 
of surrounding silicon atoms, which indeed was observed 
in EPR and ENDOR studies [24,58,59]. Theoretical calcu- 
lations which would take into account such interactions for 
neutral and positive charge states of iron are lacking. In the 
following discussion we will show that even experimentalists 
do not agree if the positive or neutral iron diffuses faster. 

The charge state of iron in n-Si and p-Si at different tem- 
peratures can be easily calculated from the position of the 
Fe? /+ level with respect to the Fermi level. The temperature 
dependencies of the fraction of ionized iron in p-Si and n-Si 
for the temperatures from 300 K to 1400 K are presented in 
Fig. 5. The calculations took into account the temperature de- 
pendence of the band gap width and of the position of the 
iron level with respect to the valence band, as shown in Fig. 4 
and discussed in Sect. 1.2.2. At room temperature, Fej is pre- 
dominantly positive (and paired with boron) in p-type silicon 
and neutral in n-type silicon. At temperatures above approxi- 
mately 100 to 200 °C the pairs dissociate, and the Fej 1 " con- 
centration increases. At sufficiently high temperatures (above 
600 K) the Fermi level moves to the midgap position, and 
iron becomes mostly neutral in both n-type and p-type sil- 
icon. The fraction of ionized iron depends on the function 
chosen to describe the temperature dependence of the iron 
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Fig. 5. The equilibrium fraction of ionized interstitial iron in p- and n-type 
sUicon for different assumed temperature dependencies of the energy level 
of interstitial iron. The simulations were made for a doping level of 
10 15 cm -3 and for an iron contamination level much lower than the boron 
(phosphorus) doping. The dotted line (see curve J in Fig. 4), the solid line 
(curve 2 in Fig. 4), and the dashed line (curve 3 in Fig. 4) correspond to 
the parabolic, steplike, and proportional temperature dependencies of the in- 
terstitial iron energy level. The calculations for p~Si also take into account 
iron-boron pairing, although the concentration of FeB pairs is not shown in 
the figure. A peak in the Fej + plot at about 400K in p-$i is due to the dis- 
sociation of FeB pairs (see Sect 2.1, 2.6 for more details). Almost all the 
ionized iron is paired with boron in p-Si at temperatures below 350 K 



level position. Parabolic and step-like functions (curves 1 and 
2 in Fig. 4, dotted and solid lines in Figs. 4, 5) result in the 
fraction of ionized iron decreasing with increasing tempera- 
ture from 10% at 1000 K to 1.5% at 1400 K in both p-Si and 
n-Si. On the contrary, the proportional model (dashed lines in 
Figs. 4, 5) predicts a gradual increase of the fraction of ion- 
ized iron from 5% at 1000 K to 10% at 1400 K. 

The predominantly neutral charge state of iron at elevated 
temperatures agrees with the conclusion of Weber [118], 
based on the theoretical model of Van Vechten [119], who 
estimated effective formation enthalpies AHf of metals in 
silicon of A// f ^ 2.4 eV for neutral interstitials and AHf& 
1.75 eV for singly positively charged interstitials. Since these 
predicted values are close to the experimentally observed for- 
mation enthalpies of Cr, Mn, Fe, Co (AJ7 f ^2.1eV) and 
Ni and Cu (A/f f ^ 1.5 eV), respectively, Weber [118] sug- 
gested that Cr, Mn, Fe, Co dissolve as predominantly neutral 
species, while Ni and Cu dissolve as positively charged ions. 
Although experimental data on the charge state of metals at 
the diffusion temperature are scarce, the available data con- 
firm the model that most of the iron is neutral at the diffusion 
temperature of 1 100 °C [120] while interstitial copper in sili- 
con is positively charged [36,120]. 

However, this does not necessarily imply that the diffu- 
sion coefficient of iron at high temperatures is determined by 
the diffusivity of neutral iron. While about 95% of iron is neu- 
tral at T > 1000K in moderately doped silicon (see Fig. 5), 
the small fraction of ionized iron may significantly change (or 
even determine) the apparent diffusivity of interstitial iron, if 
the diffusion coefficient of the ionized iron is much higher 
than that of the neutral iron. This will be discussed below. 

At low temperatures, the charge state of iron is better de- 
fined. It is assumed that iron is neutral in n-Si and in the 
depletion region of reverse-biased Schottky diodes made on 
p-Si, and is positive (and paired with boron) in p-Si. How- 
ever, generation of electron and holes by iron (see, for ex- 
ample, [9]) may make small fractions of iron ionized in the 
depletion regions of p-Si Schottky diodes, and neutral in the 
depletion regions of n-Si Schottky diodes, as pointed out by 
Heiseretal. [121], 



13.2 Diffusivity of iron: experimental data. The first expres- 
sion for the diffusivity of iron obtained by using data in a wide 
temperature range was reported by Weber [19, 122]. He made 
a fit to the data of Struthers [123], Kimerling et al. [124], 
and Shepherd et al. [125] and derived a diffusion coef- 
ficient of £>(Fei) = 1.3 x 10- 3 exp(-(0.68eV)/* B ncm 2 /s. 
Since then, more than 10 research groups reported their 
results on iron diffusivity [25,35,65,77,90,96,103,106,121, 
123, 125-139]. Some of them tried to distinguish the diffusiv- 
ity of the neutral and the positively charged iron. The diifu- 
sivity of iron at high temperatures was measured by studying 
the diffusion profile of Fej using radiotracer method [25,35, 
123, 126] or DLTS [127, 128]. Iron diffusivity at low tempera- 
tures (below 300 °C) was measured using three techniques: 
(a) temperature-dependent kinetics of association of iron- 
acceptor pairs [77,90, 125, 130-132] in p-Si; (b) precipitation 
of iron at quenched-in and extended defects [103,131,133, 
134] in p-Si and n-Si; (c) studies of iron diffusion across 
the space-charge region of a Schottky diode [121,135,136]. 
These data are plotted in Figs. 6, 7 and are discussed below. 
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Fig. 6. The diffusivity of iron in silicon at temperatures above 250 °C The 
solid line represents the fit given by (2). The fit was made using data from 
both this figure and Fig. 7 
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Fig. 7. Diffusivity of iron in silicon at temperatures below 250 °C. Solid line 
represents the fit given by (2). The fit was made using data from both this 
figure and Fig. 6. The filled diamonds and open squares are the data ob- 
tained by Koveshtkov et aJ. [135] for neutral and ionized iron, respectively 

The diffusivity data were plotted separately for the high- 
temperature and the low-temperature range to make the plots 
more legible. We want to point out that the horizontal and ver- 
tical scales of Figs. 6, 7 are identical (so that the readers may 
combine them into a single oversized graph by cutting and 
pasting) and a single line is fitted through both graphs. 

1.3.2.1 Diffusivity of iron at high and intermediate tempera- 
tures. The first four data points for iron diffusivity at high 
temperatures (7> 700 °C) were obtained by Struthers [123], 
who reported the diffusion coefficient of D = 6.2 x 10"* 3 x 
exp(— 0.87eV/JkBr)cm 2 /s. His measurements were done 



in a narrow temperature range (from 1100°C to 1265 °C) 
and featured noticeable scatter (see Fig. 6). Uskov [129] 
evaluated iron diffusivity from the kinetics of iron out- 
diffusion from the sample at 1200 °C and obtained a value 
of D = 4x 10" 6 cm 2 /s. Antbnova et al. [128] measured the 
depth profiles of Fe after annealing the samples at tempera- 
tures between 1000 °C and 1200 °C with short heat pulses 
and obtained a good agreement with the profiles simulated 
using the data of Struthers [123]. Isobe et al. [127] deter- 
mined the diffusion coefficient of iron of D = 9.5 x 10~ 4 x 
exp(— (0.65cV)/ksT)cm 2 /s in the temperature range from 
800 to 1070 °C from the depth profiles of iron distribu- 
tion measured by DLTS. Gilles obtained the diffusion co- 
efficient of D = 1.7 x 10" 6 cm 2 /s at 920°C [35,126] and 
1.4x 10~ 6 cm 2 /s at 900 °C [126] from iron out-diffusion 
experiments. 

There is only one set of data points obtained in the inter- 
mediate temperature range (between 250 °C and 400 °C) by 
Schwalbach et al. [25], who evaluated the Fa diffusion coeffi- 
cient from diffusion broadening of Mossbauer spectra of 57 Fe 
nuclei. 

7.5.2.2 Diffusivity of iron at low temperatures. The experi- 
mental data obtained in the temperature range of 0°C to 
250 °C are plotted in Fig. 7. Hie experimental conditions of 
some of the studies were such that iron is expected to be pre- 
dominantly in neutral charge state. For instance, Lebedev et 
al. [103] and Takahashi et al. [131] measured precipitation 
of Fe? in /i-type silicon using the photocapacitance method 
and EPR, respectively, and reported that there was no ap- 
preciable Fe? precipitation at temperatures less than 130 °C. 
The temperature dependence of the precipitation time con- 
stant at T> 130 °C yielded a diffusion barrier for iron of 
0.78 eV [103] and 0.80 eV [131]. However, since the density 
of precipitation sites was unknown, the pre-exponential factor 
in the expression for diffusion coefficient could not be deter- 
mined. Heiser and Mesli [121 ] used a method originally sug- 
gested by Baruch [140] to study the diffusion of Fe, in space- 
charge regions of Schottky diodes. The method consists of the 
application of a reverse bias to a Schottky diode for a long 
period of time (1 to lOOh) and monitoring the drift of the 
impurities out of the depletion region by measuring the de- 
crease of the corresponding DLTS peak amplitude. They con- 
cluded that the diffusion coefficient of neutral iron is given 
by D= 10- 2 exp(-0.84eV/fc B r)cm 2 /s (crosses in Fig. 7). 
Tavendale et al. [137] used the same technique at a single 
temperature of 25 °C and obtained D = 2.3 x 10* 16 cm 2 /s. 
Koveshinkov et al. [135] (filled diamonds in Fig. 7) moni- 
tored profiles of redistribution of iron in the depletion region 
of a Schottky diode during anneals and reported the diffusion 
barrier for Fe? of 0.56 eV. Szawelska et al. [96] monitored 
precipitation of neutral iron in w-Si and determined the dif- 
fusion activation energy of neutral iron of 0.66 ±0.01 eV. 
Finally, Lee [133] derived the activation energy of Fef diffu- 
sion of 0.69 ±0.03 eV from kinetics of precipitation of Fe? by 
EPR in the temperature range from 400 K to 470 K. 

Studies of the diffusivity of positively charged iron 
were done using predominantly p-type silicon. Heiser and 
Mesli [121] studied diffusion and drift of Fej 1 " in the deple- 
tion region of a Schottky diode and determined the diffusion 
barrier of Ftef of 0.69 eV (open triangles in Fig. 7). Their ex- 
perimental data suggested that the diffusivity of Fej** at 300 K 
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is higher than that of Fe? by a factor of 50. Koveshnikov et 
al [ 1*351 studied redistribution of iron in the space-charge re- 
irion of a Schottky diode and, in disagreement with Heiser 
et al. [121], reported that the diffusion barrier of Fe^ was 
0.92 eV (open squares in Fig. 7), i.e., higher than that of 
Fe ,j ~ and that the absolute values of diffusivity of Fe^ were 
lower than those of Fef. Boltaks et al. [134] studied precip- 
itation of Fef at dislocations and derived the diffusion co- 
efficient D(Fe?) = 6.3 x 10^exp(-(0.58eV)/* B ncm 2 /s. 
Nakashima et al. [77] determined Fe^ diffusivity from the re- 
action of association of FeB pairs and suggested the diffusion 
coefficient D(Fe^) = 3.3 x KT 1 exp(-0.81 eV/k B T)cm 2 /s. 
From similar studies, Shepherd and Turner [125] and Reiss 
et al. [141] concluded that the diffusion barrier was 0.75 eV, 
whereas Wijaranakula [132], Lagowski et al. [106], and Ober- 
meier et al. [8] obtained a good agreement with the diffusion 
barrier of 0.68 eV suggested by Weber [19]. Kimerling et 
al. [138,139] and Lemke [90] concluded from association 
kinetics of the pairs that the diffusion barrier is 0.85 eV 
and 0.83 eV, respectively. Association studies of Kaniava et 
al. [142,143] resulted in activation energies of 0.49 eV and 
0.55 eV for two different wafers. This unusiially Iqw value 
was explained by the high injection level used for the car- 
rier lifetime measurements to determine the FeB pairing 
rate. They suggested that the high density of injected nji- 
nority carriers might retard the pairing process by partial 
photo-dissociation of FeB pairs during the carrier lifetime 
measurement cycle [143]. Swanson [65] studied temperature 
dependence of decay of the quenched-in level at Ey -h 0.4 eV 
(which was unambiguously iron) and found that the activation 
energy of this process was 0.81 ±0.04 eV. 

J. 3.3 DiffUsivity of iron: analysis of data. Thus, there is asig- 
nificant scatter from 0.49 eV to 0.92 eV in diffusion barriers 
reported for interstitial iron. Additionally, there is no agree- 
ment whether ionized or neutral iron diffuses faster. While 
most of the authors that compared diffusivities of Fef and 
Fe? agree that neutral iron has a higher diffusion barrier than 
Fej* - [ 1 2 1 , 1 3 1 , 1 36, 144] , Koveshnikov et al. [ 1 35] came to the 
opposite conclusion. Heiser et al. [145] suggested that this 
contradiction can be partly explained as a consequence of 
FeB pairing, i.e., that Koveshikov et al. [135] measured the 
effective diffusion coefficient of the ionized iron, affected by 
pairing with boron. It is known that pairing may decrease the 
apparent diffusivity by orders of magnitude [146]. The impor- 
tance of a correct evaluation of pairing with shallow acceptors 
can be illustrated by the example of copper in silicon. As 
we recently showed [147], the generally accepted expression 
for copper diffusivity [36] underestimated the copper diffu- 
sion barrier by more than a factor of two, and thus the room 
temperature intrinsic diffusivity by three orders of magnitude, 
due to an incorrect accounting for the pairing of copper with 
boron. As will be discussed below, the data on the pairing 
constants of iron with boron remain ambiguous, thus hamper- 
ing accurate evaluations of the effect of pairing, as was the 
case of Cu in Si in [147]. 

It is clear that an accurate value for the diffusion barrier 
for neutral or ionized iron can be obtained only by using ex- 
perimental data in a wide temperature range. Unfortunately, 
this is difficult if not impossible since iron usually co-exists in 
both neutral and ionized charge states. The apparent diffusion 
coefficient of iron diffusing in two charge states is given by 



Dapp(Fei) == /D(Fe^) + ( 1 - /)D(Fef), where / is the frac- 
tion of ionized iron, and £>(Fen and D(F$) are the diffusion 
coefficients of ionized and neutral iron [35]. For example, if 
we assume that the diffusivity of Fe^ at a temperature T is 
much higher than the diffusivity of Fe| J (D(Fef) » D(Fef)), 
then the apparent diffusion coefficient D m (T) will be deter- 
mined by the diffusion barrier of Fef , £ ra (Fe^), despite the 
relatively low fraction, /, of the ionized iron. However, the 
pre-factor in the diffusion coefficient will be reduced by the 
value of coefficient /, which is, in turn, temperature depen- 
dent. The consequence of the temperature dependence of / is 
that neither the diffusion barrier of Fe? nor that of Fej** can be 
straightforwardly extracted from the temperature dependence 
ofD ap p(Fei). 

Despite the lack of understanding of the charge state of 
diffusing iron, the experimental data points for both Fe? and 
Fe^, all plotted in the same plot, can be fitted by a single 
straight line (Figs. 6, 7) with a slope corresponding to the dif- 
fusion barrier of £ m (Fei) = 0.67 ± 0.02 eV. The majority of 
experimental data points can thus be described by the equa- 
tion 

ZXFeO = (1.01*5) x 10- 3 exp ("^^) . cm 2 /s. (2) 

It is worth noting that this expression agrees within the error 
limits with the one previously reported by Weber [19,118], 
Z)(Fe0 = 1.3 x 10-*exp(-(0.68eV)/A B ncm 2 /s. Although 
almost all data points in Figs. 6, 7 are fitted by (2), we would 
not conclude from this result that the diffusion barrier of 
iron does not depend on its charge state. We believe that this 
agreement actually confirms the difficulties with the analysis 
of iron diffusivity due to co-existence of iron in both charge 
states in all or almost all imaginable experiments. Since we 
can assign (2) to the diffusivity of neither ionized nor neu- 
tral iron, we will call it the "effective diffusion coefficient of 
iron". Equation (2) is sufficient for most technological pur- 
poses such as prediction of intrinsic gettering behavior of 
iron. Yet, the physics of iron diffusion can not be understood 
without further experimental and theoretical studies. 

The discrepancy between the FeB pairs association ex- 
periments, which in most cases resulted in a Fef diffusion 
barrier from 0.75 eV to 0.8 eV, and the overall fit of most ex- 
perimental data to the diffusion barrier of 0.67 eV, deserves 
a separate discussion. We believe that this discrepancy is ex- 
plained by at least two reasons. 

First, the FeB association experiments could be made, at 
least partly, in a temperature range where the association rate 
of the pairs is close to the dissociation rate. In this case, the 
rate of the observed reaction is given by the sum of the as- 
sociation and dissociation rates of the pairs. An Arrhenius 
plot of the apparent reaction rate will then give an activation 
energy which lies somewhere between the barriers for associ- 
ation and dissociation of the pairs. This will be discussed in 
detail in Sect. 2.1.6. 

Another possible explanation follows from the recent ex- 
perimental studies by Takahashi et al. [148,149], who sug- 
gested that the height of the potential barrier for a jump of 
Fe[*" from the second-nearest interstitial position to the first 
nearest interstitial position to the acceptor depends on the 
chemical nature of the acceptor (0.75 ±0.01 eVfor Al, 0.72± 
0.01 eV for Ga, and 0.69 ± 0.01 eV for In), and is higher than 



the potential barrier for diffusion of Fe|*" ions from the bulk of 
the sample to the second-nearest position. This latter barrier 
was found to be 0.68 ±0.01 eV [148], in a good agreement 
with the high-temperature studies. It is interesting to note that 
the reported height of the energy barrier for the last jump 
increases with decreasing atomic number of the acceptor. Ex- 
trapolating the dependence reported in [148] to boron, one 
obtains a value of about 0.78 eV, which in fact is very close to 
the barrier for association of FeB pairs reported by a number 
of researchers. 

The diffusivity of Fej in heavily doped wafers is poorly 
investigated. We are aware of only one study, performed 
by Gilles et al. [35,126]. Using radioactive tracer tech- 
niques, they found that the diffusivity of iron at 700 °C 
is D = 4x 10" u cm 2 /s for the phosphorus doping level 
of M> = 1 x KPcm" 3 , and D = 2.4 x 10~ 7 cm 2 /s for the 
boron doping level of Afo = 8 x 10 19 cm -3 . The intrinsic dif- 
fusivity of iron at the same temperature, given by (2), is 
D = 3.4 x 10~ 7 cm 2 /s. The small decrease of the diffusion 
coefficient in p-Si was ascribed to the pairing reaction of Fej 
with B, and the pronounced decrease in n-Si was explained 
by the significant fraction of substitutional Fe in strongly 
/i-doped Si or by possible pairing reactions with P [35, 126]. 
An alternative explanation for a strong decrease of the Fe dif- 
fusivity, which follows from the discussion above, might be 
the diffusion of Fej exclusively in the neutral charge state in 
/i + -Si. This is plausible since the Fermi level remains in the 
upper half of the band gap in n + -Si even at 700 °C, and the 
positively charged fraction of iron is negligible. 



1A Solubility of iron in silicon 

The solubility of iron increases with temperature below the 
eutectic temperature, passes through a maximum, and fi- 
nally decreases above the eutectic temperature, which is 
1207 °C for the Si-rich FeSi 2 phase [150-152]. For other 
Fe-Si phases, the eutectic temperature varies from 1200 ? C 
to 1250 °C [150]. The retrograde solubility of iron above the 
eutectic temperature was reported by Trumbore [153] and 
Feichtinger [154]. The solubility of iron below the eutectic 
temperature was studied by a number of groups [35,77,83, 
84, 123, 133, 155-157]. Their experimental data are plotted in 
Fig. 8. The data published before 1980 were reviewed by We- 
ber [19,156], who suggested the following expression for the 
solubility of iron in silicon: 

5=1.8xl0 26 exp(-(2.94±0.05eV)/A:Br) l cm' 3 . . (3) 

This expression was derived using data points in a relatively 
narrow temperature range (from 950 °C to . 1200 °C) and is 
shown by a dashed line in Fig. 8. Not all data points available 
at that time were used in the linear fit since Weber [ 1 56] found 
that the solubilities of iron obtained by EPR or DLTS (open 
symbols at temperatures above 950 °C K in Fig. 8) are in most 
cases 2 to 3 times lower than the values obtained by neutron 
activation analysis (NAA). At least four possible reasons for 
that were discussed in the literature. First, it was suggested 
that part of the iron may dissolve substitutionally and may 
not be detected by techniques sensitive only to the interstitial 
iron [157]. However, as discussed in Sect 1.1, later stud- 
ies failed to detect substitutional iron in quenched samples. 
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FIg.8. Hie solubility of iron in silicon below 1200°C. The dashed line 
corresponds to the 1983 plot of Weber [19], (3), whereas the solid line is 
the new fit given by (4). The solid line was obtained by a nonlinear least- 
squares fit to solid symbols in the figure (see discussion in the text). The 
data points reported by Lee et al. [133] are represented by three different 
symbols depending on the temperature range where they were obtained, as 
well as the quenching rate used (see the graph legend and discussion in the 
text). The solubility data obtained on heavily doped silicon are not included 
in this plot and are presented in Fig. 9 



Weber [19,156] pointed out that an insufficient concentra- 
tion of iron on the silicon surface, or formation of a surface 
silicon oxide acting as a diffusion barrier for iron (see, for 
example, [158]) may decrease the dissolved concentration of 
iron below its equilibrium solubility. On the other hand, the 
concentration qf iron may exceed the equilibrium solubility 
if the equilibrium with a surface silicide is not established. 
Such a peak in the concentration of dissolved iron, exceeding 
the equilibrium solubility by up to 50%, was observed dur- 
ing the first 20-30 min of diffusion [19]. The third possible 
reason for the disagreement of experimental data are traces 
of oxygen in argon or nitrogen gases, frequently used as the 
ambient for iron diffusion. The formation of an oxide layer 
at the silicon surface, acting as a diffusion barrier, or oxida- 
tion of the iron, may result in a lower concentration of iron 
in the bulk than its equilibrium solubility [78,156,159,160]. 
This problem is usually solved by using reducing gas, a mix- 
ture of hydrogen with nitrogen or argon. Finally, the most 
important reason for disagreement in iron solubility during 
the metal diffusion is, from bur point of view, that the quench- 
ing rate used in some experiments might be insufficient to 
quench all iron interstitially and prevent the formation of 
complexes and agglomerates of iron (see discussion below in 
Sect 4, 5), the gettering of iron by oxygen precipitates or ex- 
tended defects (see [133]), or the outdiffusion of iron (see, 
for example, [121,136,161]). Colas et al. [162] even sug- 
gested that oxygen precipitates may reduce iron solubility in 
silicon. Later it was understood that they actually observed 
gettering of iron by the oxygen precipitates during cooling. 
Clearly, precipitated/agglomerated iron can not be detected 
by methods selective to Fei, and leads to ain underestima- 
tion of the total iron concentration. Summarizing the discus- 
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sion above rehable values for iron solubility can be obtained 
after anneals that are sufficiently long to establish the iron 
sihcide boundary phase (about , h)% c * nducted fa a ™ 

gas ambient, and terminated by a quench in a vertical furnace 

^i^n° f 1^ ***** ^ ^ ««««« ^ 

T> 950 C obviously did not comply with these requirements 
(for example used slow cooling rates), they resulted in under- 
estimated solubility values and were not used in the fit. These 
data points are indicated by open symbols in Fig. 8. 

Lee et al. [133] and Gilles et al. [35] reported several 
data points between 700 °C and 950 °C that lie above the ex- 
frapolation done by Weber [155,156] (open diamonds and 
triangles at T< 1 lOO K in Fig. 8). The origin of this discrep- 
ancy became clear only after the recent studies by McHuao 
et al. [83,84], who measured solubility of iron in intrinsic 
moderately P-doped (6.5 x 10 14 B/cm 3 ), and heavily boron- 

i° ped ^n 5 ^ X 10 B/cm > silicon in temperature range 
from 800 C to 1 lQO°C using NAA. In agreement with the 
previously reported data of Gilles [35], they found a signifi- 
cant increase m solubility in heavily doped substrates (Fig 9) 
McHugo et al. [83,84] showed that this solubility enhance- 
ment m p + silicon is significant only at low temperatures 
T< 1000 °C (Fig. 8). Furthermore, the dependence on the 
doping level atr< 1000°C is so strong that even a rela- 
tively low boron doping (6.5 x 10 14 B/cm 3 ) results in a sub- 
stantial increase in solubility, as shown by open circles in 
Fig. 9. Thus, the data obtained on intrinsic and doped silicon 
at temperatures below 1 000 °C should be distinguished in the 
solubility plot shown in Fig. 8. We therefore excluded from 
the following analysis several data points which were ob- 
tained at T< 1000 °C on non-intrinsic silicon (open symbols 
at T< 850 °C in Fig. 8). A linear fit made using filled symbols 
in Figs. 8, 9 yielded the following updated expression for the 
solubility of iron in intrinsic silicon in the temperature range 
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Fig.9. The solubility of iron in silicon as a function of the substrate dop- 
ing level. The solid line is given by (4) (the same line as in Fig. 8) 
and corresponds to the intrinsic solubility of iron. The doping levels are: 
1 x ltf^ cm" 3 (phosphorus), open squares; 8 x 10 19 cm' 3 (boron), open di- 
amonds; 1.5xl0 19 cm~ 3 (boron), open triangles; 6xl0 I4 cm" 3 (boron), 
open circles. The filled circles correspond to intrinsic silicon 



800 °C to 1200°C: 

, « ( 2.86±0.05eV\ _ 3 
S = (8.4!^) x 10 25 exp I ^ J , cm 3 . (4) 

It is worth noting that (4) agrees within the error limits with 
the (3) reported by Weber in 1983 [19]. A discussion of the 
formation enthalpy of iron above the eutectic temperature can 
be found in the 1983 review of Weber [19]. Those data are 
not updated in this review since, to the best of our knowledge, 
no new data on solubility of iron at silicon above the eutectic 
temperature have been reported during the last 15 years. 

The authors of this review are aware of only one theoret- 
ical calculation of the solubility of iron made by Fistul and 
Shmugurov [163]. They developed a model of potentials of 
the interaction of the impurity atom with the surrounding sil- 
icon lattice and then calculated the enthalpy of formation of 
Fef of 2.02 eV, which is reasonably close to the experimental 
data. 



2 Pairing of interstitial iron with shallow acceptors 

2.7 Iron-boron pairing 

2.J.1 EPR studies o/FeB pairs. The possibility of pairing 
of iron with boron was suggested as early as in 1957 by 
Collins and Carlson [26]. They did not, however, detect the 
pairing since their experiments were done at temperatures too 
high for a significant iron-boron pairing to occur. Evidence 
that iron-boron pairs do form was obtained from the EPR 
studies by Ludwig and Woodbury [28,45] (see Appendix A). 
They detected the paramagnetic resonance associated with 
the (lll)-oriented FeB pairs and suggested that their para- 
magnetism arises from three unpaired electrons of the 3d 7 
iron impurity in the crystal field of the negative acceptor. 
According to their model, the spin density of the ionized ac- 
ceptor is small and does not substantially contribute to the 
resonance of the pair. Although the real spin of the pair is 
5=3/2 [164,165], the symmetry of the resonance can also 
be described using the effective spin 5=1/2. The superhy- 
perfine splittings due to interstitial 57 Fe and 11 B, along with 
the very weak 29 Si ENDOR line, measurable only along the 
(111) direction, were reported by Spaeth et al. [166]. 

The resonance of the acceptor state of FeB pairs was de- 
tected by Gehlhoff and Rehse [167,168] using EPR under 
constant illumination (see Appendix A). They found that the 
spectrum is (1 1 1> axially symmetric and identified it as a FeB 
pair from the observed line splitting due to isotopes l0 B, 11 B, 
and ^Fe. 



2.7.2 Electrical, optical and recombination properties of FeB 
pairs. The identification of the energy level of (FeB)°'+ at 
Ey + (0. 10 ±0.01) eV was made by Feichtinger [157] from 
Hall effect studies and confirmed by a number of Hall ef- 
fect and DLTS studies of various research groups [80,91,92] 
Graff and Pieper [78] and Kimerling et al. [124] corroborated 
the identification of this level as that of the iron-boron pairs 
by monitoring the dissociation and association kinetics of the 
pairs. Graff and Pieper [78] also showed by using infrared 
absorption spectroscopy that the sum of the concentrations 
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of the levels at 2?v+0.38eV and Ey+QA eV remains con- 
stant during the association/dissociation reactions, and that 
the concentration of the unpaired substitutional boron de- 
creases during the association reaction and recovers during 
dissociation. The correlation between concentrations of Fej 
and FeB was confirmed by the EPR studies of Weber [169]. 
A decrease in the boron concentration and an equal increase 
in the concentration of FeB pairs was also observed in Hall 
effect studies [154], The donor nature of the FeB pairs was in- 
ferred from the comparison of the resistivity data and the Hall 
effect data before and after the pairing reaction [78,154]. 

The hole capture cross-section of the FeB pairs was re- 
ported to be CTp = 1.5 x 10" 13 cm 2 at T = 55 K [92], and 
Op = 2x 10" 14 cm 2 as measured by DLTS between 50 K 
and 65 K [80] (see Appendix B). No thermal activation of 
the capture cross-section of holes by FeB pairs was re- 
ported [80]. Zoth and Bergholz [107] reported the electron 
capture cross-section of FeB pairs at room temperature to be 
a n =4 x 10" 13 cm 2 . However, this value may be somewhat 
inaccurate since they did not take into account the possibility 
of recombination through the acceptor level of the pairs (see 
discussion below). 

Besides the well-known donor level at Ey + 0.1 eV, the 
FeB pair has also an acceptor level in the upper half of 
the band gap. Using a simple ionic model of the FeB pair, 
Feichtinger [154], Kimerling et al. [138] and Lemke [170] 
assigned the electrical level of FeB at Ey + 0.1 eV to the 
double-donor state of Fe^"*"*" , pulled up from the valence 
band into the band gap (Fig. 10) by the electrostatic inter- 
action with the nearby charge of opposite sign (Bg ). As- 
suming that the same interaction would raise the position of 
the level of interstitial iron Fej*'"*" in the negatively charged 
FeB pair by the electrostatic interaction energy of 2sj on ic = 
q 2 /(4neeor) = 0.52 eV (where r — 0.235 nm, the distance be- 
tween the closest substitutional and interstitial sites), one 
would expect [171] the acceptor state (FeB)~ to have a level 
in the upper half of the band gap at about E c - 0.25 eV 
(Fig. 10). This prediction was found to be in remarkably 
good agreement with the experimental findings of Broth- 
erton et al. [91], £ c -0.29eV, Lemke [90], E c -0.23eV, 
and Nakashima et al. [172], E c - 0.29 eV, all obtained using 
p-type samples and capacitance spectroscopy techniques with 
minority carrier injection. Gehlhoff and Rehse [167, 168] es- 
timated the energy position of the (FeB) ~ 10 level in the band 
gap at E c - (0.25 ± 0.05) eV from the observed onset of the 
spectral dependence of the photo-induced change of the EPR 
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Fig. 10. A diagram of the energy levels of the interstitial iron in the electri- 
cal field of a shallow acceptor (after [138, 154,170]). This diagram is based 
on the ionic model of iron-boron pairs 



signal. As was mentioned above, this EPR spectrum was 
identified as that of the FeB pairs from the observed line 
splitting due to the isotopes I0 B, n B and 57 Fe. The identifi- 
cation of the level at E c - 0.29 eV as (FeB)" was confirmed 
using Fourier transform spectroscopy (FITR) by Ghatnekar- 
Nilsson [173,174]. They identified four sharp lines in the 
FTIR spectrum at 6969.6, 6978.4, 7008.3, and 8291.7 cm- 1 
as the absorption by trigonal FeB pairs and evaluated from 
their experimental data the positions of two levels of FeB 
pairs at £ v + 0.1 1 eV for the (FeB)°/ + and at E c - 0.275 eV 
for the (FeB)"/° pair. 

Lemke [90] reported the capture cross-sections of (FeB) - 10 
at T = 90 K of cr n & 10" 15 cm 2 and a p ^ 3 x 10~ 15 cm 2 for 
electrons and holes, respectively. Waltz et al. [175,176] re- 
ported similar values of cr n ^ 2.5 x 10" 15 cm 2 and cr p ^ 
3 x 10" 14 cm 2 at room temperature. 

Both interstitial iron and FeB pairs are strong recombina- 
tion centers. The recombination activity of FeB pairs at low 
injection levels is about 10 times lower than that of interstitial 
iron [106,107]. The minority carrier diffusion length meas- 
ured by surface photovoltage (SPV) before (L 0 ) and after 
(Li) dissociation of the pairs is quantitatively related to the 
concentration of interstitial iron, Npe, as follows [107]: 

= 1.06 x 10 16 x (Lf - L~ 2 ) (cm" 3 ) , (5) 

where Lq and L\ are measured in microns. It is some- 
what unexpected from the Shockley-Read-Hall (SRH) statis- 
tics [177,178] that the level of FeB pairs at £ v + 0.1eV, 
being so shallow, is such a strong recombination center. This 
controversy was explained [91] by the suggestion that re- 
combination, particularly at low injection levels, takes place 
via the acceptor state at Eq - 0.29 eV which is closer to the 
midgap than the donor level of the FeB pairs. This model was 
later confirmed by a number of studies of the temperature de- 
pendence of minority carrier lifetimes in iron-contaminated 
silicon [142, 179] and the dependence of the effective lifetime 
on minority carrier injection level [142,175,176]. A com- 
plete model of recombination via the FeB pairs, valid for 
both high and low injection levels can be developed only 
taking into account both the donor and the acceptor levels 
of FeB [180,181]. As it can be shown from the SRH statis- 
tics (see [182]), the recombination activity of shallow levels 
increases as compared to that of deep levels with increas- 
ing injection level rj = An//?o, where po is the equilibrium 
density of holes in the sample, and An is the density of in- 
jected electrons. A number of groups [142,143,181,183,184] 
showed that the recombination activity of FeB pairs actually 
becomes higher than that of Fej if the injection level is high. 
Thus, it is important to remember that (5) is valid only for 
low injection levels as used in the SPV technique. Further- 
more, the coefficient 1 .06 x 1 0 16 in (5) is valid for typical 5 to 
10 Q cm silicon and may vary depending on the Fermi energy 
level, i.e., on resistivity of the wafers [185] 

Only the trigonal (111) symmetry has been detected 
for FeB pairs, whereas orthorhombic (100) configurations 
have also been reported for the other donor-acceptor pairs 
(Sect. 2.2, 2.3). It is important to emphasize that die differ- 
ence between the trigonal and the orthorhombic configura- 
tions of the iron-acceptor pair is not only in.the symmetry, but 
also in the distance between the iron and the acceptor. This 
distance is r = 0.235 nm (the nearest-neighbor distance be- 
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tween Fej and B s ) for the trigonal symmetry of the pair and 
r = 0.272 nm (the second nearest-neighbor distance between 
Fej and Bs) for the orthorhombic symmetry of the pair. In 
some cases, we will also have to mention the second nearest 
position along the (1 1 1) direction. This actually corresponds 
to the third nearest-neighbor distance between Fe ; and B s . 
r = 0.450 nm. 

A number of researchers tried to detect the orthorhom- 
bic configuration of FeB pairs using capacitance spectroscopy 
methods. Dobaczewski et al. [86,186] applied the high- 
resolution Laplace-transform DLTS technique to the capaci- 
tance signal of FeB pairs and found a well-defined single peak 
in the Laplace spectrum, which corresponds to a monoexpo- 
nential emission from the FeB level at Ey +0. 1 eV. They also 
found that another level with a hole emission activation en- 
ergy of 0.074 eV appeared in the spectrum after dissociation 
of FeB pairs by minority carrier injection at 150 K. The level 
was found to be unstable at room temperature and annealed 
out after several hours. As with the pairs of levels associ- 
ated with FeAl, Fein, and FeGa complexes (see Sect. 2.2, 
2.3, 2.4, 2.6 below), it was suggested that this level is an un- 
stable donor level of the orthorhombic configuration of FeB 
pairs [86]. The position of this second level of FeB at 
0.07 eV is in agreement with the general trend specific to 
iron-acceptor pairs that the level of the orthorhombic config- 
uration of the pair is closer to the valence band edge than that 
of the trigonal configuration. 

A whole set of deep levels associated with different con- 
figurations of FeB pairs has been reported by Nakashima 
et al. [172,187-189], who studied dissociation of FeB 
complexes under low level minority carrier injection CI to 
100mA/cm 2 ) at T— 150 K. The idea of their experiment 
was that the minority carrier injection will assist iron to 
make a jump from the interstitial position nearest to boron 
to one of the next neighbor positions. If the temperature is 
low enough to prevent thermally activated diffusion of iron, 
these metastable configurations can be 'frozen in" and will 
become detectable by electrical measurement techniques. In- 
deed, they found four unstable electron traps at Ec — 0.43 eV, 
E c - 0.46 eV, E c - 0.52 eV, E c - 0.54 eV, and two hole traps 
at £ v 4-0.53eV and £ v +0.48eV. These levels annealed 
within lOmin at temperatures ranging from 200 K to 240 K. 
The annealing kinetics [187, 189] revealed activation energies 
for each of the traps varying from 0.65 eV to 0.75 eV with 
pre-exponential factor close to 10 12 -10 13 s" 1 as would be ex- 
pected for a single-jump process (i.e., the lattice vibration 
frequency). The activation energy was explained [1 87, 1 8 8] as 
the barrier to atomic motion of Fej 1 " ion from one configura- 
tion to another. The concentration of the FeB acceptor level 
at Ec - 0.29 eV, which was the only level in the upper half of 
the band gap before the minority carrier injection, decreased 
in concentration after injection by a value equal to concentra- 
tion of the emerging metastable levels. It was suggested [ 1 88] 
that the level at E c - 0.43 eV is the acceptor level of the ( 1 OO) 
orthorhombic configuration of FeB pairs. 



2.1.3 Equilibrium binding energy of FeB pairs. The forma- 
tion and dissociation reactions of FeB pairs are determined 
by three energies schematically presented in Fig. 1 1 : the equi- 
librium binding energy £ b » the diffusion barrier of interstitial 
iron Ejn, and the dissociation energy of the pairs iidiss- 




Fig. 11. An energy diagram of the FeB interaction in silicon (after Kimer- 
ling et a!. [138]). E d ^ is the dissociation energy of the pair, E b is the 
equilibrium binding energy (a difference in equilibrium energy m the paired 
and unpaired states), and E m is the diffusion barrier for iron. Note that 
£ diss < £b + £ra. and that there is no straightforward manner to calculate 
one of the three parameters from the other two 



The equilibrium binding energy E b can be obtained from 
the analysis of the temperature dependence of the equilibrium 
fraction of Fei, paired with boron, as it follows from the law 
of mass-action: 

^expf-^V (6) 

where N\ is the density of interstitial sites in the silicon lattice 
(5 x 10 22 cm" 3 ), and Z is the number of possible orienta- 
tions of the pair with the same symmetry around one ac- 
ceptor atom (4 for the tetrahedral symmetry). Note that (6) 
defines the equilibrium ratio of concentrations of FeB pairs, 
//(FeB), to the concentration of the ionized iron, 7/(Fei + ), 
rather than to the total concentration of unpaired iron. The 
value of W(Fei + ) is usually smaller than the total concentra- 
tion of the unpaired iron, the ratio of the two being dependent 
on the doping level of the sample and the temperature. It is 
also important to keep in mind that the equilibrium binding 
energy determines the paired fraction of iron at equilibrium* 
i.e., after infinitely long equilibration time. The kinetics of 
the association/dissociation reactions are determined by is m 
and £dj SS . 

The temperature dependence of the equilibrium fraction 
of paired iron was studied by Kimerling et al. [138,139], 
Lemke [90,170], Reiss et al. [141], and Wunstel et al. [80] . 
Hie reported binding energies vary from 0.45 eV [139] and 
0.53 eV [80] to 0.6eV [141,170] and 0.65 eV [90,138]. Wi- 
jaranakula [132] calculated a binding energy of 0.58 eV 
from the experimental data reported earlier by Nakashima 
et al. [77]. A scatter in binding energies reported by differ- 
ent groups can be partly ascribed to interaction of iron with 
other impurities, unintentionally introduced in silicon during 
studies. A careful re^arnination of the pairing data using 
state-of-the-art silicon material and modem cleanroom facil- 
ities is required to obtain the binding energy with a better 
accuracy. Our own experience suggests that the true value 
of the binding energy lies in the range from 0.58 eV, as re- 
ported in [77,132,170], to 0.65 eV, as reported in [90,138]. 
That these values are close to the electrostatic binding en- 
ergy fjonic = 0.52 eV of a donor-acceptor pair is often use<I 
as an argument in favor of the model of ionic type binding of 
the FeB pair. However, a number of theoretical calculations 
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discussed in Sect 2.6 strongly suggest that the binding be- 
tween Fej and shallow acceptors, in particular boron, has not 
only ionic character, but also has either a covalent or a strain 
component. 

2.1 A Association kinetics ofFcB pairs. The association ki- 
netics of FeB pairs can be described [125, 190] using the theo- 
ries of diffusion-limited precipitation/trapping, developed by 
Ham [191], and of trapping of mobile donors by immobile ac- 
ceptors, developed by Reiss et al. [146]. In the latter theory, 
a mobile ion is considered captured when the electrostatic at- 
traction potential between the donor and the acceptor, 7(r), 
exceeds the average thermal energy, k B T y i.e., k B T < V(Rc). 
Since in most cases Rc is as large as several nm, the interac- 
tion potential can be approximated by Coulombic attraction 
even for the mixed type of binding in the pair, and Rq is then 
obtained by solving the following equation: 

where p is the free hole density. For the doping levels 
Na < 10 17 cm" 3 and under extrinsic conditions, the screened 
Coulomb potential in (7) can be approximated by a simple 
Coulomb potential (i.e., the exponential term in (7) can be 
neglected), and the capture radius Rc can be obtained explic- 
itly as Rc = q 2 /(4n8£oksT). Inserting this expression into 
the formula for the diffusion-limited precipitation time con- 
stant (Ham's law, [191]), which for die case of trapping of 
iron ions by boron has the form: 

ta SS = (4^D(Fe i )^ A /?cr I , (8) 

one can easily obtain [106, 124,138,190] 

ee 0 k B T _ 566.7x7 
855 s 2 tf A 0(Fei) ~ Z>(Fei) x N A ' 

where Na is the concentration of shallow acceptors (in our 
case, boron) in cm -3 , T is the temperature in Kelvin, and 
D(Ft{) is the iron diffusivity in cm 2 /s. Equations (8) and (9) 
assume that the interstitial iron concentration is much lower 
than that of the boron. Hieslmair et al. [190] showed that 
Ham's law [191] can be used to describe the association ki- 
netics and obtained a capture radius of 5 to 6 nm at room 
temperature, in a good agreement with the theoretical pre- 
diction of approximately 4.5 nm from (7). Equation (9) was 
extensively used to determine the diffusion barrier of ion- 
ized iron at low temperatures. These studies were discussed 
in Sect. 1.3.2.2 and yielded strongly scattered values for the 
diffusion barrier from 0.58 eV to 0.8 1 eV. 

2.1.5 Dissociation of FeB pairs. The dissociation kinetics of 
FeB are determined by the potential barrier E^ss for a jump of 
theFei ion away from the first closest-neighbor position to the 
boron atom. The dissociation time constant is then given by: 

-K^)- <"» 

where v is the attempt frequency. Although it is quite 
straightforward to determine the dissociation barrier E&ss 



from the temperature dependence of dissociation reaction, 
we are aware of only two values reported in the litera- 
ture: £ diss = 1.17eV (v=1.8xl0 10 s- 1 ) [154], and 1.2eV 
(v = 5 x 10 10 s~ l ), reported by the same author in [157]. In 
the latter article, Feuchtinger studied "migration" of the FeB 
pairs, which clearly consists of dissociation of a pair, diffu- 
sion of iron to the next ionized boron, and re-association of 
the pair. The kinetics of the reaction is thus determined by 
the dissociation reaction. We want to emphasize that we did 
not find any indication that the process of re-association of 
iron-boron pairs (see Sect 2.1.6) was properly taken into ac- 
count during the dissociation studies in [154,157]. Thus, it 
can not be excluded that the dissociation energy may need to 
be re-evaluated. 

Ostapenko and Bell [192,193] reported the dissociation 
of FeB pairs by using ultrasound. At 75 °C, about 25% of 
donor-acceptor pairs could be dissociated by ultrasonic treat- 
ment. It was suggested that the ultrasound vibrations can 
force a partial separation of Fe\ and B. This reduces the FeB 
pair binding energy and accelerates the process of complete 
pair dissociation. 

The dissociation reaction of FeB pairs is also enhanced 
by minority carrier injection. Kimerling and Benton [138] re- 
ported that the dissociation reaction could be driven to greater 
than 95% completion by the injection of minority carriers at 
temperatures above 77 K. The dissociation rate of the pairs 
was reported to increase linearly with the injection current 
and was weakly temperature dependent with the activation 
energy of about 0.1 eV [139]. Minority carriers can be in- 
jected either by forward biasing p~n junctions (and even 
aluminum Schottky diodes), or by shining white light on the 
wafer. Light-stimulated dissociation of FeB pairs was first re- 
ported by Graff et al. [78] and is now routinely used in SPV 
systems [106]. The dissociation was carried out by focus- 
ing the light of a 50 W halogen lamp onto a 0.5-cm-diameter 
spot on the surface of the wafer [78]. The minority car- 
rier injection-stimulated enhancement of dissociation of the 
pairs is explained as a recombination-enhanced defect reac- 
tion [194]. Examples of recombination-enhanced diffusion in 
silicon are enhancement of diffusion of interstitial aluminum 
in silicon [195] and the decrease of the barrier for migra- 
tion of iron from the first-neighbor position in FeAl, Fein, 
and FeB pairs to the second one [196-198]. Reaction rates 
that have been increased by 6 orders of magnitude using light 
were observed in other systems such as GaAs [199-201], 
GaP [202,203], and GaAsi_^P x [204]. Weeks et al. [205] 
suggested that the energy released during the capture of a mi- 
nority carrier at the defect and the recombination with trapped 
majority carrier is entirely converted to vibrations of the de- 
fect, reducing the apparent activation energy of the purely 
thermal annealing behavior or the diffusion of the defect. The 
attempt frequency of the reaction would correspond to the 
electron-hole recombination rate. The observed reduction of 
the activation eneigy of dissociation of about 0.8 eV would 
fit to a model of converting the electron capture energy of 
0.77 eV on Fe^ to defect vibrational modes. Further theor- 
etical studies of nonradiative capture and multiphonon emis- 
sion [206,207] (see also [208] for a review of non-radiative 
transitions in semiconductors) enabled Sumi [209,210] to de- 
velop a model of recombination-enhanced defect reactions 
that takes into account energy dissipation to the surrounding 
lattice. 
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2.1.6 FeB pairing: analysis. An important factor which is 
frequently neglected in analysis of reactions of FeB pairing 
and dissociation is that the rate of the observed reactions is 
determined by both the association rate = and the dis- 
sociation rate r d iss = t^. In most studies the association or 
dissociation energies of the pairs were calculated from the 
temperature dependence of the reaction rate assuming that the 
reaction is either purely dissociation or purely association. 
This assumption is in most cases incorrect and may lead to 
substantial errors. This point can be illustrated by the analy- 
sis of a simple differential equation, describing the reaction 
kinetics. If we neglect the precipitation of iron, assume that 
most of the interstitial iron is ionized (W(FeO = //(Fe^)) and 
that the concentration of iron is much less than the concentra- 
tion of boron, then the iron— boron pairing reaction is given by 
the following differential equation: 

d7V(FeB) 

— = rass x (AXFei) - N(FcB)) - r diss x //(FeB) , 

(11) 

where W(FeB) is the concentration of iron-boron pairs, and 
N(Ft\) is the concentration of ionized interstitial iron. The 
solution of this equation is given in the general case by 

tf(FeB) = r °* MFeO-l r ~ N(Fe i )-Afo(FteB)l 

Toss + r^ss J Tags + f*diss J 

x exp(-(r ass +r diss )0 , (12) 

where // 0 (FeB) is the initial concentration of FeB pairs at the 
beginning of the measurement (t = 0). It is important that the 
reaction rate is always given by the sum of association and 
dissociation rates. Thus, the temperature dependence of the 
reaction rate is not necessarily determined by only the dis- 
sociation energy or only the diffusion barrier. The apparent 
activation energy of the reaction as reported in the literature 
may actually come out to be between the true dissociation and 
diffusion barriers. In this case, the apparent energy would de- 
pend on the doping level of the samples and the temperature 
range used in the studies. A complete description of the re- 
action kinetics should include the temperature dependence of 
the ionized fraction of iron [132, 170,21 1,212] and the precip- 
itation (agglomeration) rate of iron. 

2J2 Iron-aluminum pairing 

2.2,1 EPR studies of FeAl pairs. Unlike iron-boron pairs, 
which are stable only in (1 1 1) configuration, iron-aluminum 
pairs co-exist in two configurations. One of the first EPR stud- 
ies of FeAl pairs was that of van Kooten et al. [213], who 
detected resonances associated with the trigonal and ortho- 
rhombic configurations of (FeAl) 0 pairs. Both spectra (see 
Appendix A) showed a sixfold splitting caused by the hyper- 
fine interactions of iron with one Al nucleus. An additional 
twofold splitting appeared in the resonance of the orthorhom- 
bic (FeAl) 0 pair in the samples doped with iron enriched to 
90% in the magnetic isotope 57 Fe. This clearly indicated the 
involvement of one iron and one aluminum atom in the para- 
magnetic center. The association and annealing behavior of 
the center is qualitatively similar to FeB pairs. After keeping 
the samples for one night at room temperature, the resonance 



of Fei disappeared and resonances of (FeAl) 0 pairs increased 
in intensity by about a factor of 3. Annealing at T = 90 °C 
for 1 h was found to be sufficient to dissociate the pairs [213]. 
Gehlhoff et al. [214] reported that there are actually two sets 
of EPR lines associated with the orthorhombic (FeAl) 0 pairs. 
They showed that these lines are caused by the spin transi- 
tions within an isolated quartet ground state with a zero-field 
splitting very large compared to the microwave quantum, and 
suggested to describe the lines using a single Hamiltonian 
with a spin S = 3/2 [214,215] (see Appendix A). 

The initial EPR studies of van Kooten et al. [213] and 
Gehlhoff et al. [214] discussed above revealed only the res- 
onance of the neutral charge state of (1 1 1) FeAl pairs. EPR 
resonances of negative, (FeAl)", and positive, (FeAl) + , pairs 
with a trigonal symmetry have been found only recently 
by Gehlhoff and Irmscher [168,215,216] (see Appendix A). 
Whereas the spectrum of a neutral pair (FeAl) 0 could be 
easily observed in all samples containing Fe and Al, the ap- 
pearance and intensity of the other two spectra were much 
more dependent on the sample preparation and illumination. 
The resonances associated with (FeAl) 0 and (FeAl) + were 
identified from the observed isotope splitting [166,215]. For 
the (FeAl)" pair such a straightforward proof was not pos- 
sible because line broadening prevented the resolution of the 
hyperfine structures. It should also be noticed that Irmscher et 
al. [215] suggested that the trigonal (FeAl)-related spectrum 
reported earlier by Greulich- Weber et al. [217,218] and ten- 
tatively identified as (FeAl) 0 complex with the Fe; + ion at the 
second nearest position from the Al atom along the (1 1 1) di- 
rection, is identical to the resonance, interpreted as (FeAl) + 
in [2 1 5] . Therefore, we listed both resonances in the same row 
in the table in Appendix A. 



2.2.2 Electrical levels of FeAl pairs. The electrical levels 
of the (111) and (100) configurations of FeAl pairs were 
detected in electrical measurements even before they were 
observed by EPR. Feichtinger [154] and later Lemke [170] 
and Nakashima et al. [219] found from Hall effect and DLTS 
measurements that the concentrations of the electrically ac- 
tive Fej*" and Al§ decreases within several hours after dif- 
fusion and quench of Fei h* Al-doped silicon. Additionally, 
a level of FeAl pairs at Ey -4- (0.20±0.01) eV appears in the 
band gap. The second energy level of FeAl pairs at Ey -f- 
(0.13±0.01)eV was reported in [78,80,144,170,219,220], 
It was shown that the pairing of interstitial iron with alu- 
minum results in an increase of the concentration of the levels 
of FeAl and in a decrease of the concentration of the level 
of Fe;, such that the sum of concentrations of all three lev- 
els remained constant [78]. Optical measurements using the 
photocapacitance method were in good agreement with the 
optical and thermal ionization energies of the level Ey -h 
0.20 eV [81]. Hole capture cross-sections were reported to 
be cr p = 2 x 10~ 14 cm 1 [220] and a p = 4 x 10" 15 cm 2 [80] 
for the level at Ey + 0.20eV, and cr p = 6 x 10 -14 cm 2 [220] 
for the level at £V + 0.13eV. The equilibrium binding en- 
ergy of the FeAl center with £v + 0.20eV was reported to 
be 0.70 eV [170]. An analysis of the dependence of the hole 
emission rate on electric field in the depletion region of 
a Schottky diode revealed no significant shift in the position i 
of the FeAl DLTS peaks [144] (as compared to the shift of j- 
14 K, predicted by the Pool-Frenkel law for an acceptor level ^{ 
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in p-type silicon [221]), confirming that both configurations 
of FeAl pairs are donors. 

Theoretical calculations for FeAl pairs are sparse. Weih- 
rich and Overhof [222] predicted from ab initio calculations 
that the pair with orthorhombic symmetry should be stable in 
the negative and in the neutral charge states. Trigonal pairs 
were found stable in negative, neutral, and positive charge 
states and in two structural configurations: with Fej on the 
nearest and with Fej on the second-nearest interstitial pos- 
ition. 



2.2.3 Metastability <?/FeAl pairs. The metastability of FeAl 
pairs was discovered in 1985, when Chantre et al. [144] as- 
sociated the FeAl energy levels at £v + 0.20eV and Ey + 
0.13 eV with the trigonal and orthorhombic configurations 
of FeAl pairs, respectively. They also showed that the pairs 
can rearrange from trigonal to orthorhombic configuration 
and back under the influence of an applied bias or light. 
This was explained as a single diffusion jump of Fei from 
a nearest neighbor position to substitutional aluminum (that 
corresponds to the trigonal symmetry of the pair) to a second- 
nearest neighbor position (that corresponds to the ortho- 
rhombic symmetry of the pair). They found that the relative 
amplitude of the two DLTS peaks of the FeAl pair depends 
on the cooling conditions of the sample before the DLTS 
scan, as shown in Fig. 12. Both DLTS peaks were present 
when the sample was cooled with a reverse bias. Cooling 
with a zero bias was found to eliminate the signal of the trap 
at J?y + 0.13eV, whereas the signal from the trap at Ey + 
0.20 eV grew correspondingly. Hie kinetics of transitions 
between the two defect states revealed that the rate of transi- 
tion from the metastable (Ey + 0.13 eV) to the stable (Ey + 
0.20eV) state is given by 1 x 10 u exp(-(Q.5Q6V)/k B T) 
(s _1 ), whereas the rate of the reverse reaction is given by 
5 x 10 12 exp(-(0.64eV)/* B r) (s _l ). 
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Fig. 12. The DLTS spectra of FeAl pairs in silicon. The amplitudes of the 
DLTS peak are associated with the stable and metastable states of the FeAl 
pair, and depend on the cooling conditions of the sample (after Chantre et 
al. [144]). The DLTS rate window is 1000 s^ 1 



Chantre and Bois [144] obtained a good quantitative 
agreement with the experimental data assuming that the in- 
teraction between iron and aluminum is fully described by 
the electrostatic attraction between them. Gehlhoffet al. [223] 
established unambiguous correlation between the deep lev- 
els detected by DLTS and EPR resonances. They showed that 
the metastable state can be populated not only by cooling 
down the sample with bias, but also by illuminating it dur- 
ing cooldown with above-band-gap light. This was explained 
as recombination-enhanced conversion from one to the other 
configuration of the pair. 

More details of the metastable properties of iron-acceptor 
pairs can be found in Sect. 2.6. 

2.3 Iron—gallium pairing 

The EPR resonance parameters of the neutral (FeGa) 0 pair 
with a trigonal symmetry were reported by Ludwig and 
Woodbury [28,56]. The orthorhombic configuration of the 
FeGa pairs was detected by Gehlhoff et al. [214]. They de- 
tected two sets of EPR lines and proved that they are as- 
sociated with the spin transitions within an isolated quartet 
ground state of the same center with a zero-field splitting very 
large compared to the microwave energy. The resonance was 
correlated with FeGa pairs from the well-resolved hyperfine 
structures [214]. Even the hyperfine splitting caused by the 
isotopes ^Ga and 71 Ga could be distinguished both for the 
orthorhombic and the trigonal pair [214,224]. 

DLTS studies of electrical levels associated with iron in 
gallium-doped silicon initially revealed only one FeGa level 
at Ey + 0.24 eV [80, 170] (see Appendix B). The equilibrium 
binding energy of this pair was reported to be 0.47 eV [170]. 
Optical measurements using the photocapacitance method re- 
vealed a level at £v + 0.24eV [81], in good agreement with 
the DLTS measurements. The hole capture cross-section of 
the FeGa center with a Ey + 0.24 eV level was reported to be 
cr p = 2 x 10" 14 cm 2 by Wtinstel et al. [80] and a p = 3 x 10~ 15 
(at 120 K) by Lemke[170]. 

The bistability of iron-gallium pairs, similar to that of 
FeAl pairs, and the second energy level of FeGa pairs, was 
reported by Chantre and Kimerling [225]. In addition to 
the previously reported level of FeGa pairs at Ey + (0.23 ± 
0,01) eV [80, 170], they reported a metastable level at Ey + 
0.14eV [225]. The elastic energy barriers for the transform- 
ation from the metastable to the stable position and back were 
found to be about 0.5 eV and 0.6 eV, respectively. 

2.4 Iron-indium pairing 

2.4.1 EPR identification of Fein pairs. The EPR signal of 
Fein pairs in silicon with a (100) symmetry was reported in 
1962 by Ludwig and Woodbury [28] (see Appendix A). Later 
experiments by Gehlhoff et al. [226,227] using isotope split- 
ting confirmed the assignment of this spectrum to Fein pairs. 
In contrast to FeB, FeAl, and FeGa pairs, the stable config- 
uration of Fein pairs was found to be orthorhombic. EPR 
identification of the metastable trigonal configuration of Fein 
pairs was made by Omling [228,229] (see Appendix A). The 
metastable trigonal Fein signal was shown to appear in the 
EPR spectrum of /?-type samples only if they were illumi- 
nated with photons of hv > 0.5 eV. Later, one more reson- 
ance associated with the trigonal Fein complex was found by 
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Gehlhoff et al. [226,230,231]. From the temperature depen- 
dence of the intensity of the EPR signals, they concluded that 
the spectrum originally detected by Ludwig and Woodbury 
originated from an excited state, while the newly discovered 
resonance is consistent with the ground state of the pair. 

An EPR resonance tentatively assigned to the positive 
charge state of the orthorhombic Fein center was found by 
Gehlhoff et al. [227,232] (see Appendix A). The negative 
charge states of the Fein complexes have not yet been found 
by EPR, but have already been detected by FTIR, as discussed 
below. 

2.4.2 Electrical and optical properties of Fein pairs. The 
electrical levels of the stable orthorhombic and metastable 
trigonal configurations of Fein pairs were detected using 
DLTS by Chantre and Kimerling [225], who reported Fein- 
related levels at Ey + 0. 1 5 eV and E w + 0.27 eV. The elastic 
energy barrier for the transformation from metastable state 
to the stable configuration and back were found to be about 
0.5 eV and 0.6 eV, respectively. 

Tidlund et al. [233-235] studied excitation spectra of Fein 
pairs using FTIR and found a Fein-related set of two lines 
at about 6100 era" 1 760 meV) and a considerably weaker 
set of 6 lines at about 6600 cm- 1 (^820meV). They as- 
sociated these sets of lines with the orthorhombic and trig- 
onal configurations of Fein pairs, and estimated the pos- 
itions of the (FeIn)~ /0 acceptor levels to be E c — 0.39 eV 
and E c - 0.32 eV. 

Reislohner [236] applied perturbed angular correlation 
spectroscopy (PAC, see [237]) to Fein complexes and found 
three different configurations of Fein pairs. One of the config- 
urations corresponded to the (1 1 1) orientation of the pair; the 
second one was tentatively assigned to the (100) orientation. 
The third configuration did not show any definite orienta- 
tion dependence and was assigned to a complex with a large 
number of equivalent orientations, such as a pair in a (31 1} 
orientation. All pairs annealed out at T> 200 °C. The upper 
limit of the dissociation barrier of the complex was estimated 
to be £ D (FeIn) < 1 .4 eV [236]. 

In contrast to other Fe-acceptor pairs, calculations for 
the Fein pairs based on simulations of the potential energy 
profiles using Lennard-Jones potentials and electrostatic in- 
teractions, suggested that the second-nearest to In interstitial 
position of Fej in the (1 1 1) direction is more stable than the 
nearest-neighbor position for the neutral and positive charge 
state of the pair [238]. Since the interaction with the 29 Si lig- 
ands could not be resolved by EPR measurements [229] and 
since both interstitial sites are consistent with the observed 
trigonal symmetry, it is not possible to distinguish these two 
arrangements by EPR. However, Gehlhoff et al. [224] argued 
that since the trigonal crystal fields estimated from the meas- 
ured g- values [239] are approximately the same for FeAl and 
Fein, the assumption that Fe ion can be in the second intersti- 
tial position in (1 1 1) direction seems to be very improbable. 

2.5 Iron-thallium pairing 

Thallium is one of the possible /?-type dopants in silicon. 
However, to the best of our knowledge, there are no data 
on electrical levels or binding energies of FeTl pairs, pri- 
marily due to the scarcity of thallium-doped silicon sam- 
ples. The data on iron-thallium complexes are limited to 



a few photoluminescence (PL) studies [240-242]. The pos- 
ition of the PL lines associated with FeT] was reported to 
be around 1.086 eV, measured at liquid helium temperature. 
The Si:Tl luminescence at 1 .086 eV was found to vanish com- 
pletely above 12 K, being replaced by a similar set of lines at 
1 .049 eV. This led the authors of [242,243] to suggest that the 
complex can exist in two different configurations. The trig- 
onal structure of the Si:Tl luminescent center, suggested by 
Sauer and Weber [244], was later confirmed by means of de- 
tailed stress and Zeeman studies by Watkins et al. [245] and 
Conzelmann et al. [246]. They also showed that the high- 
temperature center shows a reduction in symmetry, possibly 
to an orthorhombic configuration with (100) symmetry. 

The intensity of the PL lines in Si:Tl spectrum was greatly 
enhanced if the samples were heated to at least 800 °C and 
then quenched. Even a short anneal for several seconds was 
sufficient to activate the centers [241,242,247]. Schlesinger 
and McGill [248] reported that the intensity of the PL lines 
in Si:Tl samples increases with increasing Fej concentration, 
which can be considered an indication that iron is a part of 
this complex. 

Despite the apparent agreement on the identification of the 
Fell PL signal in the above-mentioned references, we would 
exercise caution with respect to these data since the identifi- 
cation of the PL spectra initially attributed to the other, better 
investigated iron-acceptor pairs remains highly controversial 
and even contradictory, as it will be discussed in Sect. 6. 



2.6 General trends in iron - shallow acceptor pairing 

2.6.1 General picture of pairing and dissociation of iron- 
acceptor pairs. Although there is a wide disagreement be- 
tween different research groups about the binding energies 
of the iron-acceptor pairs, the qualitative picture of their as- 
sociation and dissociation kinetics is well established. After 
a fast quench, iron is predominantly interstitial and is strongly 
supersaturated since its equilibrium solubility at room tem- 
perature is vanishingly small. Shallow acceptors are effect- 
ive trapping sites for iron around room temperature in p-Si 
due to electrostatic attraction between the negatively charged 
acceptors and the positively charged iron ions. Thus, iron- 
acceptor pairs are formed within several hours to several days 
at room temperature depending on the average distance be- 
tween the ionized iron and the acceptor atoms. Annealing of 
the sample at temperatures above 200 °C dissociates the pairs. 
As observed from EPR and DLTS measurements, complex- 
ing of iron in p-type silicon exhibits three phases. The tem- 
peratures of the transition between phases are similar for all 
p-type dopants [81,91,132,148,149,174,211,249]. In phase 
I, at temperatures below 80 °C to 100 °C, the concentration of 
interstitial iron either decreases as the pairs are formed, or re- 
mains at a very low level. In phase n, at temperatures between 
100 °C and 200 °C, a recovery of the Fe* species and a re- 
duction in the concentration of the Fe~acceptor pairs is ob- 
served [250] (see Fig. 5). The temperature at which the max- 
imum concentration of ionized interstitial iron is achieved 
depends on the doping level [249]. Hie pairing and disso- 
ciation processes in phases I and II are reversible. In phase 
III, at temperatures above 200 °C, acceptors can not pro- 
vide stable traps for iron and iron diffuses towards more 
stable sinks such as surface defects or dislocations where it 
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precipitates/agglomerates [52]. The concentrations of both 
Fej and Fe-acceptor pairs decrease as the isochronal anneal- 
ing temperature increases and do not recover after the end of 
the anneal [142], 

2.5.2 Type of binding in the pairs. The role of electrostatic 
and covalent bonding in iron-acceptor pairs is one of the most 
intensely discussed iron-related fundamental topics. There 
is extensive experimental evidence that the electrostatic at- 
traction between Fej 1- and negatively charged substitutional 
shallow acceptors is a driving force for the formation of iron- 
acceptor pairs. This was confirmed by the comparison of the 
pairing rate of iron in p-Si and in n-Si partly compensated 
with boron. No pairing was reportedly observed in bulk of 
n-Si even after a year [90,103,138]. Similarly, no pairing 
was observed in the space-charge region of a reverse-biased 
Schottky diode on p-Si where interstitial iron is neutral even 
after 15 days [90]. These observations are frequently used as 
an argument in support of the ionic type of binding in the 
pairs. Ludwig and Woodbury [28] applied the ionic model 
to explain the EPR resonances of iron-acceptor pairs. They 
suggested that the resonance of the metal-acceptor pair is as- 
sociated with the 3d shell of the metal in the electrostatic field 
due to the negatively charged acceptor. Although this model 
was quite successful in explaining EPR spectra of the pairs, 
the assumption that the complex stabilizes as an ionic pair 
was questioned by Assali, Leite et al. [251-257]. In theor- 
etical calculations of the electron structure of iron-acceptor 
pairs made using a 25Si + B s Fei cluster and based on the 
multiple-scattering Xa formalism within the framework of 
the Watson-sphere-terminated molecular-cluster model, they 
concluded that the covalent bonding effects play a basic role 
in determining the physical properties of iron-acceptor pairs 
and include, besides Fe and B, the silicon atoms nearest to 
the pair. In agreement with Assali et al [252-255], Sugi- 
moto and Seki [258] showed from ab initio molecular orbital 
model calculations that iron loses 0.37 electrons during for- 
mation of FeB pairs. However, this charge is transferred to 
the surrounding silicon atoms rather than to the substitutional 
acceptor, which implies that there is no direct covalent bond 
between iron and boron. Instead, Fe is weakly bound to the 
silicon atoms next to boron. In contrast to the above men- 
tioned studies, recent ab initio total-energy computations of 
Overhof and Weihrich [259,260] did not reveal any covalent 
bonding between Fet and B s . In the framework of their model 
Weihrich and Overhof [260,261] also calculated parameters 
of the hyperfine interactions for iron-acceptor pairs, and ob- 
tained a good agreement with the experimental data. It is not 
yet clear what the reason is for this major disagreement in the 
theoretical calculations discussed above. 

2.63 Two configurations of the pairs and their bistability, 
EPR studies discussed in Sects. 2.1.1, 2.2. 1, 2.3, 2.4 revealed 
that Fe-acceptor pairs have two structural configurations, 
with trigonal and orthorhombic symmetry. Each of these con- 
figurations may exist in negative, neutral, and positive charge 
states. This results in up to 6 different EPR resonances (which 
may each consist of two groups of lines corresponding to 
upper and lower Kramers doublets) and up to 4 deep levels 
for each Fe-acceptor pair. Interstitial iron occupies the first 
nearest-neighbor interstitial sites adjacent to acceptors along 
the (111) direction in the pairs with trigonal symmetry, and 



the second nearest-neighbor interstitial sites adjacent to ac- 
ceptors along the (100) direction in pairs with orthorhombic 
symmetry. Both orientations of the pairs were observed at 
room temperature in Al, Ga, and In-doped samples, whereas 
the orthorhombic configuration of the FeB pair was only ten- 
tatively identified by capacitance spectroscopy methods after 
minority carrier injection at 7= 150 K [86,172,187-189]. 
The experimental observation that FeB, FeAl, and FeGa pairs 
are stable in the trigonal configuration, whereas Fein prefers 
orthorhombic configuration is in agreement with the calcu- 
lations of Kimerling [238], who explained the instability of 
(111) configuration for Fein as a result of competition be- 
tween repulsive (elastic) and attractive (electrostatic) energy 
terms, and with the ab initio molecular orbital method calcu- 
lations of Sugimoto and Seki [258]. 

The co-existence of two possible structures of the pairs 
results in the bistability of the Fe-acceptor complexes. The 
metastability and bistability of defects is a well-known phe- 
nomenon in semiconductor physics [262-265]. The two 
configurations of the Fe-acceptor pairs can be viewed as 
two energy wells separated in a configuration space by 
a potential barrier (Fig. 13). While this model can be ap- 
plied to all iron-acceptor pairs, the energies are different 
for different acceptors. The experimentally determined po- 
tential barriers separating the stable and metastable con- 
figurations of FeB, FeAl, FeGa, and Fein pairs in both 
neutral and positive charge states (see Fig. 13), and the bar- 
riers calculated using the electrostatic model of the pairs, 
are presented in Table 1. The configuration with the low- 
est total energy will have the highest population at equi- 
librium. Relative populations of the states with the (111) 
and (100) symmetries are given by A^((lll»/iV((100)) = 
exp(A£(0)/fcB7*) for the neutral charge state of the pairs, 
and N((\ 1 1))///(<100)) = exp(AE(+)/k B T) for the positive 
charge state. Since AE(0) ^ A£(+), the relative population 
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Fig. 13. A configuration-coordinate diagram of the bistable iron-acceptor 
pairs in silicon. The FeB, FeAl, and FeGa pairs have a lower sum of the 
electronic and elastic energies in {111) configuration than in (100) configu- 
ration (as shown in the figure), whereas Fein pairs have a lower equilibrium 
energy in the (100) configuration than in the (111) configuration. The nota- 
tions in the figure are explained in the text and in the caption to Table I 
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Tablel. The parameters of the bistable states of complexes of iron with shallow acceptors, calculated using the ionic model and ^^^^S^^ 
The notation follows that of the configuration coordinate diagram Fig. 13 and is as follows: E b is the binding energy of the pairs, E(a) and E{b) are the 
hole emission activation energy in two states of the defect ((III) and (100) orientation, respectively); A£(0) and A£(+) are the differences ; in ^ jbnum 
energies of the (100) and (111) configurations of the iron-acceptor pair in neutral fFe-A)° and positive (Fe-A) + charge states respective^ £ b (0) and 
£ b (+) are the potential barriers for the jump of a positively charged Fe^ (total charge state of the pair Fej*"A^ is neutral) and d°uUy ^mvely^ed 
Fef+ (total charge state of the pair Fe?+Ar is positive) ions from the metastable to the stable position. The value of £ b (0) = 0 62eV predicted from^ the 
ionic model, was calculated as £ b (0) * - AE(+)/2, where E m is the migration barrier of Fe+. E h (+) can not be predicted from the ionic model since 
the diffusion coefficient of doubly positively charged iron is unknown 



Pair 



Ionic model 
FeB 



FeAl 



FeGa 
Fein 



£ b /eV 



£<a)/eV«IU)) 



£(*)/eV ((100)) 



A£(0)/eV A£(+)/eV £ b (+)/eV 



£b(0)/eV 
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0.04 [225] 
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of the stable and metastable configurations can be altered 
by controlling the charge state of the pairs. This is done in 
DLTS experiments by changing the bias applied to the sam- 
ple during cool down. The relative concentration of the pairs 
in each of the two configurations, and the values of AE(0) 
and AE(+) can be evaluated from the height of DLTS peaks. 
Assuming a purely electrostatic interaction between Fej and 
acceptors, one easily calculates the characteristic energies 
AE(0) and A£(+) of the pair [144] as 0.071 and 0.142eV, 
respectively. These values are in relatively good agreement 
with the parameters of FeAl and FeGa pairs (Table 1). In con- 
trast, the simple ionic model fails completely in describing 
the energetics of the Fein and FeB systems [225]. Apparently, 
in these cases either the strain associated with the oversized 
Ins or the undersized Bs, or covalent bonding, play a signifi- 
cant role in determining the defect properties. 

2.6 A Energy levels of iron-acceptor pairs. The donor en- 
ergy levels of iron-acceptor pairs were discussed in Sect. 2.1- 
2.3 and are summarized in Table 1 and in Appendix B. Two 
trends are observable. First, the (111) configuration of the 
pair always has a deeper level in the band gap than the (100) 
configuration. Hie second trend is that the levels become 
deeper with the increasing tetrahedrally coordinated covalent 
radius of the acceptors (88 pm for B, 126 pm for Al, 126 pm 
for Ga, and 144 pm for In, [266]), the ionic radius of the 
acceptors (35 pm for B, 39 pm for Al, 47 pm for Ga, and 
62 pm for In, [267]), and the increasing ionization energies 
of the shallow acceptors (B (44.3 meV) to Al (68.5 meV), Ga 
(72.7 meV), In (156 meV), [81]). The dependence of the en- 
ergy levels of the pairs on the acceptor ionization energy was 
fitted by a straight line by Lemke [170], who also pointed 
out that a similar dependence was observed for Ma-acceptor 
and Cr-acceptor pairs. We found, however, that a much bet- 
ter linear fit can be achieved if the donor energy levels of 
iron-acceptor pairs are plotted versus covalent radius of the 
acceptors, as presented in Fig. 14. Recent theoretical calcula- 
tions, which took into account the lattice relaxation due to the 



difference in size of shallow acceptors and that of the silicon 
atoms, enabled Zhao et al. [268] to obtain a fair agreement of 
the calculated energy levels with the experimental data. The 
importance of lattice strain, particularly with respect to the 
pairs of Fej with the largest of acceptors, In, was also pointed 
out earlier by Benton [263]. She suggested that the stability of 
the Fein complex in the configuration with iron in the second 
nearest position can be explained as a competition between 
the attractive electrostatic forces and the repulsive strain 
effects. 

Theoretical calculations of the energy position of the 
acceptor levels of Fe^-acceptor pairs [257], using multiple- 
scattering Xa formalism predicted the values of Eq — 0.20, 
£ c -0.14, and E c -0.10eV for FeB, FeAl, and FeGa pairs 
in the (100) symmetry, respectively. The values of £c - 0.32 
and £ c -0.31eV were predicted for the (100) configura- 



u 
m 

*c5 

Q. 

(D 

.c 



i 

Q> 
C 
LU 



U.3U 


i i - 


i i r 




0.25 




Al y 








• / 




0.20 








0.15 


















B / 






0.10 








0.05 


1 i 


— i 1 * 


* 



80 90 100 110 120 130 140 150 
covalent radius of the acceptor (pm) 

Flg.14. The dependence of the hole ionization energy on the covalent radii 
of shallow acceptor impurities. Th* filled and open circles correspond to the 
upper and die lower donor levels of the iron-acceptor pairs, respectively 
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tions of the FeAl and FeGa pairs, respectively [257], These 
values agree quite well with the experimental data presented 
in Appendix B. 

3 Pairing of interstitial iron with shallow donors 

Phosphorus is the most common n-type dopant, and its reac- 
tions with iron, if they exist, are very important in order to un- 
derstand the behavior of iron in n-type silicon and phosphorus 
diffusion gettering of iron [269]. The literature data on iron- 
phosphorus pairing are contradictory and inconclusive. Over- 
hof and Weihrich [259] inferred from ab initio total energy 
calculations that the trigonal FejP s pair is stable in n-type 
samples in the negatively charged state Fe?~Pj. They found 
that the pairs should be diamagnetic and therefore EPR inac- 
tive. However, there is no convincing experimental evidence 
for the existence of such pairs. Takahashi et al. [270,271] re- 
ported that the concentration of the interstitial iron detected in 
n-Si after diffusion and slow cool decreased with increasing 
phosphorus doping level. This was interpreted as formation of 
stable FeP pairs with the binding energy of 0.9 eV [270,27 1]. 
Conversely, Lee et al. [272] reported that isochronal anneal- 
ing for 30 min at temperatures of up to 500 °C did not change 
the intensity of the EPR line of the P donor even after Fe? dis- 
appeared completely at 200 °C. Therefore, they argued that 
Fef does not have a high probability of forming pairs with 
phosphorus. Hohne et al. [273] noted that in samples contain- 
ing a sufficiently high concentration of phosphorus or arsenic 
(more than roughly 10 15 cm -3 ), the Fef EPR line exhibited 
a satellite structure of four lines. The structure of the satellites 
was found to be isotropic, and their intensities were found to 
increase with increasing P concentration. Doping of the same 
samples with 57 Fe isotope resulted in the splitting of both 
the central FeP line and the satellites. A spin-spin interaction 
between iron and phosphorus (or arsenic) could explain the 
positions of the satellites. However, the spin-spin interaction 
was not as large as usually observed for close partners. To ex- 
plain this phenomenon, Hohne et al. [273] suggested that they 
observed a local random or almost random distribution of Fe? 
with respect to shallow donors, which implies that close pairs 
of Fe? and shallow donors do not form [273]. 

It is appropriate to emphasize that no electrostatic attrac- 
tion, which is the main driving force for the formation of 
iron-acceptor pairs, can arise between iron and phosphorus. 
This is because phosphorus is positively charged whereas 
iron is either neutral or positive. Studies of phosphorus diffu- 
sion gettering also do not provide indications of the existence 
of stable FeB pairs. While the formation of phosphorus- 
impurity pairs is considered one of the physical processes 
involved in phosphorus diffusion gettering, evidence that such 
a pairing does exist was obtained primarily for gold [274- 
277], which is known to form both donor and acceptor levels 
in the band gap. 

4 Complexes of iron with other impurities 

4.1 FeAu complexes 

The first reports on the formation of FeAu pairs were made 
by Feichtinger [157] and Lee et al. [272] in 1979. Fe- 
ichtinger [157] concluded from his EPR studies that the 



resonances of both Fef atoms and FeB complexes are strongly 
suppressed in Au-diffused samples. Lee et al. [272] found 
a new EPR spectrum with a trigonal symmetry which ap- 
peared after a 30-min anneal at 200 °C of Fe and Au co- 
diffused samples. The same resonance was also detected in 
a sample left at room temperature for a month. The spectrum 
was an order of magnitude stronger in p-type samples than in 
n-type samples and appeared only after cross-contamination 
of samples with iron and gold. The detailed EPR studies of 
Kleinhenz et al. [278] and Sieverts et al. [279] confirmed the 
existence of hyperfine interactions with the Au nucleus and 
revealed an additional twofold splitting of each line due to the 
hyperfine interaction with ^Fe nuclei [278]. The FeAu cen- 
ter could be annealed out at T = 300 °C [250,272,278]. Since 
gold is known to dissolve in silicon substitutionally, the FeAu 
pair is thought to consist of a substitutional gold atom and an 
interstitial iron atom. 

Studies of n- and p-type silicon samples co-diffused with 
Fe and Au using conductivity and Hall effect techniques 
revealed a level at about Ey + (0.4 to 0.5) eV, which was 
tentatively assigned to the FeAu complex [157], The same 
level at Ey + 0.434 eV was observed in n-Si by Brotherton et 
al. [280], who used photocurrent DLTS, a modification of the 
minority-carrier spectroscopy, described in [281]. This level 
could also be detected in p-Si. Measurements using a com- 
bination of DLTS and EPR [280], and DLTS alone [282], 
confirmed that the changes in the concentration of FeAu pairs 
during anneals anticorrelate with the concentration of gold. 
DLTS studies confirmed that the intensity of the FeAu DLTS 
peak could be substantially increased by a 30-min anneal at 
T = 200 °C to 250 °C, whereas anneals at T> 350 °C disso- 
ciated the pairs [272,280]. 

Photo-EPR of n-Si:Fe:Au proved the existence of an ac- 
ceptor level of these pairs in the upper half of the band 
gap [283]. DLTS investigations of the same sample revealed 
a level of FeAu pairs at £ c -0.33eV. The annealing ki- 
netics of the DLTS line matched that of the EPR reson- 
ance. Lemke [284] also observed the FeAu related level at 
E c - 0.35 eV (a n = 3 x 10~ 15 cm 2 at T = 160 K). Brother- 
ton [280] reported very close values of Ec — 0.354 eV, cr n = 
2 x 10 -15 cm 2 , measured at T = 80 K to 160 K. A Fe-related 
level with the same activation energy in n-Si was reported 
by Kimerling [124], although no connection to gold was 
established. 

Hie first theoretical model of the FeAu center, suggested 
by Kleinhenz et al. [278] and Sieverts et al. [279], was based 
on the ionic model of Ludwig and Woodbury [28], According 
to this model, the observed EPR signal comes from a mag- 
netic coupling between the angular momenta of gold and iron. 
It was speculated [278] that a much higher electronegativity 
of gold compared to that of iron makes the transfer of one 
electron from the iron to the gold likely, making the center 
Fei + Aus~. Strong antiferromagnetic coupling between Au<T 
(5= 1) and Fe/* (5 = 3/2) should result in a total effect- 
ive spin of S = 1/2, in agreement with experimental data. 
However, a theoretical analysis by Assali et al. [285,286] 
suggested that the FeAu pairs are formed by a covalent mech- 
anism that includes the silicon neighbors in addition to Aus 
and Fej-. According to their model, the S == 1 /2 spin of the 
center is related to a molecular orbital spread out over the 
entire cluster rather than being derived from the interactions 
between two localized magnetic centers. 
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The experimentally observed energy levels are in a very 
good agreement with the energy levels at Ev + 0.45 eV and 
E c -0.35eV, predicted by Overhof and Weihrich [287] for 
the trigonal (FeAu) 0/+ and (FeAu)~ /0 , respectively, from 
self-consistent ab initio total energy calculations. The same 
authors [259,287] predicted a covalent binding energy of the 
(FeAu) 0/+ pair of about 1 .2 eV, and about 1 .0 eV for the nega- 
tively charged pair. Theoretical calculations for a hypothetical 
FeAu pair with orthorhombic symmetry resulted in a bind- 
ing energy lower than that of the trigonal pair by 0.3 eV. This 
difference in energy would result in about 1000 times lower 
equilibrium concentration of the orthorhombic pairs at room 
temperature. 

The dissociation energy (ZT d iss in Fig. 1 1) of FeAu pairs 
was determined in the temperature range T = 350 °C to 
400 °C by Weber and Wiehl [283], who found that w = 
4.7 x \O- 3 exp{(0J3tV)/k E T). 

Experimental evaluations of the equilibrium binding en- 
ergy of FeAu pairs proved to depend drastically on the model 
used for the pairing reaction. Lemke [284] studied equilib- 
rium concentrations of unpaired Fej, Aus, and complexes 
FeAu in the temperature range between 250 °C and 400 °C in 
/i-type Si and obtained the dependence 



#(AuFe) x Nj 
//(Au s ) x //(FeO 



:exp 



/ 0.75eV \ 
\ k B T ) * 
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His calculations were based on a formulation of the mass- 
action law (13) that does not account for the different charge 
states of iron and gold, i.e., assumes that the pairing reaction 
is charge-state independent. Brotherton et al [288] argued 
that since the driving force for formation of FeAu pairs is the 
electrostatic attraction between Fe^ and AuJ , the product of 
the total concentrations of gold and iron #(Aus) x Af(Fe0 in 
the mass-action law (13) should be substituted by the concen- 
trations of the ionized impurities N(Au s ~(T)) x N(Fei + (T)\ 
where T is the temperature. This treatment resulted in a bind- 
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ing energy of 1.22±0.02eV. We want to emphasize that the 
difference in the binding energies reported by Lemke [284] 
(0.75 eV) and Brotherton et al. [288] (1.22 eV) stems from 
the different interpretation of the experimental data rather 
than from a difference in the data points themselves, as 
can be seen in Fig. 15. In this figure, we plotted the ratio 
//(AuFe)/MAu t otai)^(Fe t o ta i) from the work of Lemke [284] 
and the same ratio calculated using the data of Brotherton et 
al. [288]. There is a reasonably good agreement between both 
authors, with the exception of three data points obtained by 
Brotherton et al. [288] at low temperatures, where the kinetics 
of the reaction are so slow that the equilibrium concentration 
of FeAu pairs might not be achieved. Since the acceptor level 
of Au (Ec — 0.55 eV) lies above the donor level of Fe (£ v -f 
0.385 eV), for any position of Fermi level either Fe or Au or 
both of them will be neutral. It is easy to calculate that only 
a small fraction of either gold or iron (less than 20% at 650 K, 
and less than 1% at 400 K, [288]) will be ionized at any tem- 
perature. Low fractions of Fej*" and AuJ will not negate the 
pairing reaction between iron and gold. As long as some per- 
centage of gold and iron is ionized, the reaction will go to 
equilibrium, albeit slowly. 

A question remains whether a pairing reaction between 
neutral iron and neutral gold is viable. Such a reaction is only 
possible if the binding in the pairs is covalent, as suggested 
in [259,287]. Although the capture radius of neutral Fe by 
a trapping site (for example, by gold) (about 0.235 nm, the 
nearest distance between interstitial and substitutional atoms 
in Si) is about 20 times less than that of the Fe^ (about 5 nm, 
see [190]), the density of neutral iron and gold are much 
higher than that of the ionized ones. Thus, reactions between 
both charged and uncharged impurities may eventually need 
to be taken into account in a complete model. 

4.2 FeAg, FePd, FePt, and FeZn complexes 

Iron-silver complexes are very poorly investigated. To the 
best of our knowledge, there are no experimental data on 
electrical levels or on the microscopic structure of FeAg 
pairs, although theoretical calculations predict that these pairs 
should be stable. Overhof et al. [287] predicted two electri- 
cal levels for FeAg pairs in the band gap, at Zsv + 0.19eV 
for (FeAg)°/ + , and E c - 0.27 eV for (FeAg)"/ 0 from ab ini- 
tio calculations. They also found that the positive, neutral, 
and negative charge state of the pair should have a spin 
S = 0, S = 1/2, and S = 1 , respectively. In the following pa- 
per, Overhof et al. [259] concluded that the type of binding in 
FeAg pairs is predominantly covalent. 

The electronic structures of FePd and FePt complexes 
were studied by EPR by van Oosten et al. [289] (see Ap- 
pendix A). After diffusion of iron into platinum-doped sili- 
con, the EPR signal of platinum was no longer observed, but 
a strong signal of trigonal FePt centers appeared. The FePt 
resonance exhibited large hyperfine satellites due to the 19S Pt 
isotopes. After co-diffusion of Pt and 57 Fe, a well-resolved 
hyperfine splitting was observed. A similar trigonal spec- 
trum with hyperfine splitting due to ^Fe was observed in 
palladium-diffused samples. Alves et al. [290] have recently 
reported that an intensive 777 meV luminescence line was 
observed in n-type samples known to have some Fe con- 
tamination and speculated that this line- may be due to FePt 
complexes. 
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Czaputa [291,292] studied DLTS spectra of Si samples 
after Pd and Fe diffusion and established a correlation of 
Fe and Pd concentrations with the level at E c — 0.32 eV, 
which he attributed to FejPd s pairs. He found that anneals 
at T = 110°C could convert the palladium-related peak at 
E c - 0.22 eV into the peak at E c - 0.32 eV, whereas anneals 
at around 300 °C stimulated the reverse reaction. The activa- 
tion energy for dissociation of this FePd complex was found 
to be about 1.5 eV. Co-diffusion of Fe and Pd in p-type sili- 
con resulted in a new DLTS peak at Ey -f 0.52 eV, which has 
similar formation and annealing kinetics as the FePd complex 
detected in n-Sl No detailed studies of this defect were re- 
ported; it was only suggested that this defect may be due to 
a different charge state of the same defect [292]. Hohne and 
Juda [293] reported that the co-diffusion of iron and platinum 
in silicon suppressed the formation of platinum complexes 
consisting of 6 platinum atoms [294], although they could not 
provide any explanation for this phenomenon, 

FeZn complexes were reported, to the best of our know- 
ledge, only in [295] by Lemke, who found two new energy 
levels with £ v + P.17eV and £ c -0.47eV in DLTS spec- 
tra obtained after diffusion of iron in n- and /Mype silicon 
co-doped with zinc. He also reported that the level at Ec — 
0.47 eV did not dissociate after annealing at 200 °C for 5 min 
and suggested that it is formed by doubly negatively charged 
zinc and singly positively charged iron: Fe + Zn 2 ~. 



4.3 FeS complexes 

Three EPR resonances attributed to FeS pairs were re- 
ported by Ludwig [296] in 1965. The resonances observed in 
phosphorus-doped, undoped, and boron-doped silicon start- 
ing material were labeled FeS®, FeS(II), and FeS(III), re- 
spectively. From the dependence of their intensity on the 
doping level of the samples, Ludwig [296] suggested that 
FeS(I), FeS(II), and FeS(m) are associated with the negative, 
neutral, and positive charged states of the same FeS complex. 
The presence of Fe was proven by doping with 57 Fe, whereas 
the participation of S was inferred from the fact that the FeS 
pairs did not appear in crystals not contaminated with sul- 
fur. Later, Schirmer [297] reported one more FeS-related EPR 
signal, FeS(IV). They found that after a fast quench, only the 
EPR resonances of Fe?, P°, S + , and (S-S) + were observed. 
After several days of storage at room temperature, the Fej* in- 
tensity decreased by about a factor of 100 and the FeS (IV) 
lines appeared in the spectrum. Also the (S-S) + intensity 
decreased, pointing out that the signal FeS(IV) may be due 
to S— Fe— S complexes. The center could be dissociated by 
above-band gap light. To the best of our knowledge, there are 
no data on the energy levels of these centers in the literature. 



4.4 Interaction of iron with oxygen and carbon 

Czochralski-grown (CZ) silicon contains oxygen and car- 
bon in very high concentrations (between approximately 
cm 3 and 10 18 cm" 3 ). Although the understanding 
of the interaction of iron with carbon and oxygen is of major 
importance for the characterization of iron in silicon, there are 
only a few experimental studies on this subject which are con- 
tradictory. For example, there is evidence that iron influences 



oxygen precipitation. However, there is a contradiction on 
whether iron enhances or retards precipitation. The studies of 
Hackl et al. [298] showed that iron prevents the nucleation of 
oxygen precipitates even if the iron concentration is as low as 
10' 1 cm" 3 , whereas Zhang et al. [299], Shen et al. (300] and 
Jablonski et al. [301] found that the precipitation of oxygen 
is significantly enhanced in the presence of iron, which could 
be a result of stress relaxation around the oxide by the forma- 
tion of FeSi2- Miyazaki et al. [4] showed that grown-in iron 
enhances the formation of oxygen-induced stacking faults 
(OSF) and, following the conclusions of Hourai et al. [302] 
and Sadarnitsu et al. [303], who observed similar effects on 
wafers contaminated with iron from the surface, suggested 
that iron precipitates serve as nucleation sites for OSF. 

You et al. [304] reported an electron trap level at Eq — 
0.36 eV due to the iron-vacancy-oxygen complex (FeOV) in 
irradiated samples contaminated with iron. The identification 
of this level was based on the fact that the anneals of the 
samples at 80 °C resulted in an increase in the concentration 
of the level at Ec — 0.36 eV, accompanied by an equal de- 
crease in the concentration of the A -center at £c — 0.17eV. 
The latter level is known to be associated with V-0 com- 
plex [305—307]. The total concentration of the A center and 
the center at Ec — 0.36 eV remained constant within 1 0% dur- 
ing the whole anneal, thus indicating that the first center is 
converted into the second one [304|. No reduction in the A- 
center concentration was detected during similar anneals in 
samples containing no interstitial iron, thus confirming that 
iron is a part of the Ec - 0.36 eV center. A decrease in the 
concentration of the A -center during iron precipitation was 
also observed in the EPR studies of Kustov et al. |308]. Am- 
merlaan [57], with a reference to a private communication 
with J.W. Corbett, reported ^-values of the EPR signal of 
FeOV center (see Appendix A). 

Mchedlidze and Matsumoto |3(WJ studied the electrically 
detected magnetic resonance (EDMR) of iron in Czochralski- 
grown (CZ) silicon. They found a new signal with an ampli- 
tude correlating with the iron contamination level in the sam- 
ples. Since the signal was only observed in iron-contaminated 
oxygen-rich CZ-Si, and was not detected in low-oxygen- 
content float zone (FZ) Si, it was proposed that the signal 
originates from complexes containing both iron and oxygen, 
for example, iron-decorated oxygen precipitates. This sig- 
nal could not be detected in the EPR measurement mode, 
probably due to the lower sensitivity of EPR as compared 
to EDMR. They also reported that the trigonal EDMR sig- 
nal (see Appendix A), detected after the contamination of CZ 
samples with iron, changed to that of the orthorhombic sym- 
metry (see Appendix A) after annealing the samples at 450 °C 
for 60 h. 

Wtinstel and Wagner [79| reported a defect level at Ey + 
0.33 eV, which appeared only after a slow cool (air cool or 
furnace cool) of the samples after the iron diffusion. They 
assigned this level to an iron-oxygen complex because it ap- 
peared in higher concentrations in CZ-grown -silicon with 
a higher oxygen concentration than in FZ-grown silicon. 

These few experimental results comprise the current in- 
formation on the interaction of Fej with oxygen, and is by 
no means complete or self-consistent. The understanding of 
toe interaction of iron with carbon in silicon is even poorer. 
The only article we are aware of that establishes a corre- 
lation between iron and carbon in silicon is the paper of 
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Kustov et al. [308], who studied the compiexing of iron in 
electron-irradiated silicon with different carbon concentra- 
tions. They found that the concentration of interstitial iron de- 
creases after irradiation and reaches a constant value of about 
1/3 of the starting concentration for a low carbon density 
Nq = 5x 10 I6 cm~ 3 , and 2/3 of the starting iron concentra- 
tion for high carbon density Nq = 4 x 10 17 cm" 3 . However, 
they did not infer compiexing of iron with carbon from these 
experiments. The experimental findings were explained by 
the formation of Fe-I and Fe— V complexes in concentra- 
tions depending on the reaction rate of the vacancies (V) and 
the silicon self-interstitials (/) with carbon [308]. 

4.5 Interaction of iron and its complexes with hydrogen 

The ability of hydrogen to passivate electrically active de- 
fects is well known (see, for example, a review of Pearton et 
al. [310]). Hydrogen passivation of grown-in defects is com- 
monly used in solar cell technology to increase the minority 
carrier diffusion length (see, for example, [31 1,312]). Passi- 
vation is usually implemented by a relatively low-temperature 
hydrogen plasma treatment. Substantial concentrations of hy- 
drogen are also introduced into the near-surface region of sili- 
con during routine chemical etching (see, for example, [313- 
316]). Recent experimental data showed that interaction of 
hydrogen with impurities is more complicated than just pas- 
sivation. Complexes of transition metals with hydrogen often 
introduce new levels in the band gap (see, for example, re- 
cently published data on Pt— H complexes [3 1 5], Ag— H com- 
plexes [317], and Co-H complexes [318]). 

Despite the importance of the problem of interaction of 
iron with hydrogen and passivation of iron-related defects, 
the data in the literature are inconclusive. Tavendale and 
Pearton [137,319] studied the influence of hydrogen plasma 
treatment on Fe-related deep levels and reported that the lev- 
els Ey + 0.32 eV and E v + 0.39 eV (note that the latter level 
is not the level of Ffej, since interstitial iron had in [137,319] 
a DLTS peak at higher temperatures with the energy E\ + 
0.40 eV), which appeared in DLTS spectra after diffusion of 
Fe and were thus suggested to be complexes of Fe with an 
unknown structure. These levels were neutralized by hydro- 
gen to the depth of the hydrogen penetration. They also found 
that the level of interstitial iron was not affected by the plasma 
treatment. The latter result has recently been confirmed by 
Weber [320]. Data of Kouketsu et al. [321,322] indicate that 
the FeB peak at E\ + 0. 1 eV is effectively neutralized by hy- 
drogen. Yakimov et al. [323] suggested that the introduction 
of hydrogen during the wet chemical etching stimulates the 
dissociation of FeB pairs in the near-surface layer and in- 
creases the local concentration of unpaired interstitial iron. 
However, the mechanism of such a dissociation process re- 
mains unclear. Kaniewska et al. [324] reported a change in the 
shape of the DLTS spectra of n-type Fe-contaminated silicon 
samples with dislocations and stacking faults after hydro- 
genation. However, the complex nature of the spectra made it 
impossible to identify the passivated centers. 

Several authors suggested that the levels at Ey + 0.23 eV, 
£ v + 0.38eV [321,322], and £ v + 0.31 eV [325] appear in 
the band gap after hydrogenation of iron-contaminated sam- 
ples and argued that they are levels of iron-hydrogen com- 
plexes. It was shown that these levels are unstable at el- 
evated temperatures and dissociate after anneals at 175 °C 



for 30 min [325]. Unfortunately, the experiments described 
in [321,322,325] lacked a systematic analysis of the inter- 
action of hydrogen with the impurity levels, such as that 
recently developed in [313,317]. This approach includes the 
comparison of DLTS spectra of cleaved (no hydrogen) and 
etched (hydrogen in the near-surface region) samples, and 
a comparison of the depth profiles of the traps with the pre- 
dicted distribution of hydrogen. This enables one to even 
determine the number of hydrogen atoms in a complex. To 
the best of our knowledge, applications of this method to iron 
have not yet been reported in the literature. 



5 Complexes that include several iron atoms 

In the text above we discussed complexes of one iron atom 
with other impurities. In this section we will give a review of 
another class of iron-related defects, which consist of several 
(two to four) iron atoms. These defects were found by EPR 
and can be considered as the first stage of iron precipitation. 

The first report on such defects was published by Muller et 
al. [29], who studied complexes of iron in electron-irradiated 
silicon. Since irradiation of iron-doped silicon created a num- 
ber of new iron-related EPR resonances, it was concluded 
that these defects contain iron and irradiation-induced de- 
fects, most probably, vacancies. We want to emphasize that 
whereas the number of iron atoms was established from the 
hyperfine splitting, the number of vacancies in these com- 
plexes is somewhat speculative [29]. Spectrum NL20 (see 
Appendix A) consisted of two iron atoms and was tentatively 
identified as (FeiFejV)" complex. Spectrum NL21 was only 
observed in n-type silicon and was assigned to a (FejFej V2) 
complex [29]. Spectrum NL25 also contained two iron atoms 
and was tentatively identified as another charge state of the 
defect NL20, namely (FeiFeiV) + . However, one year later 
Ammerlaan [57] changed the designation of the NL25 spec- 
trum to belong to the FeiFei complex. 

The EPR spectrum of FejFej pairs (NL24) was identi- 
fied by Muller et al. [29] and further studied by van Kooten 
et al. [326] and Gehlhoff et al. [223]. It was found [326] 
that this complex can be observed in quenched non-irradiated 
samples, thus indicating that the complex does not contain 
irradiation-induced defects (vacancies or self-interstitials). 
The maximum amplitude of the EPR resonances was after the 
quench and disappeared after a storage at room temperature 
for 12 h, or after an anneal at 75 °C for 1 h [326]. The dis- 
appearance of this center at or just above room temperature 
indicates that there is a strong preference for the further ag- 
gregation of iron at this center. Since the EPR spectrum could 
only be fitted using a spin-Hamiltonian with a half-integer 
spin, the possibility of two neutral Fep atoms was ruled out 
Since the spectrum was observed in high resistivity n- and 
p-type silicon and in low-resistivity p-Si, it was suggested 
that the resonance corresponds to the positively charged com- 
plex (FejFei) + . This fits well with the assumption that the 
real spin of the system is S = 5/2 [326]. Unfortunately, the 
authors of [326] did not succeed in finding a unique least- 
squares fit to their experimental data for S = 5/2 and had to 
restrict their analysis to an effective spin S = 1/2. 

Gehlhoff et al. [223,224] found a resonance of a de- 
fect with a monoclinic symmetry, that was identified as an- 
other configuration of the FeiFei center. Although the authors 
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of [223] used an effective spin 5= 1/2 to determine the 
g-tensor, they reported that 5 = 5/2 corresponded to the more 
.realistic value of g & 2.05; yet they did not report a complete 
analysis with this spin. The FejFej center was found to be un- 
stable at room temperature. According to the data of [223], 
the disappearance of the EPR resonance of the FeiFej pairs 
could be enhanced by the illumination of the sample with 
above-band-gap light. 

The next intermediate state of the agglomeration of iron, 
which was observed by EPR, is the (FcOa cluster. Muller et 
al. [29] reported that the trigonal EPR spectrum NL22 (see 
Appendix A), observed in both irradiated and non-irradiated 
samples, could be modeled only with an unusually high 
effective spin 5 = 4. To explain this spin, they suggested 
a model of four ferromagnetically coupled spins 5 = 1 of neu- 
tral interstitial iron atoms Fef. The complex is thus neutral. 
Studies of the dependence of the concentration of this cen- 
ter on the temperature of isochronal anneals revealed that the 
concentration of the complex reaches its maximum at tem- 
peratures around 140 °C. If the sample is annealed at a higher 
temperature (above 150 °C), the intensity of NL22 decreases. 
This can be explained as a growth of (Fej)4 clusters to larger 
agglomerates of iron. Finally, the remaining NL22 centers 
dissociate at T> 250 °C [29]. Note that iron is not the first 
metal in silicon that was found to form clusters of four atoms. 
As early as in 1959, Ludwig and Woodbury [327] reported 
a cluster of four manganese atoms in silicon and observed 
a hyperfine interaction of the Mn atoms in the cluster. Un- 
fortunately, such an interaction could not be observed in [29] 
due to the approximately ten times smaller nuclear g-f actor of 
iron compared to manganese. 

Ezhevskii et al. [57,328,329], Gehlhoff et al. [223,224], 
and Irmscher et al. [215] reported complexes with monoclinic 
and orthorhombic symmetries consisting of two iron atoms 
and a shallow acceptor such as boron or aluminum (see Ap- 
pendix A). The centers exhibited quadropole splitting due to 
the n B isotope, confirming the involvement of one boron 
atom. When the iron impurity was enriched with the ^Fe 
isotope, an additional splitting was observed. The low sym- 
metry of the complexes and the high value of the real spin 
5 = 5/2 (note that the complete analysis of the symmetry of 
these spectra was made for the effective spin 5=1/2, see Ap- 
pendix A) suggested that the complexes consist of two iron 
atoms [328]. The reaction producing the complexes of FeFeB 
took place during storage at room temperature, and the inten- 
sities of the spectra were found to be a function of the storage 
time [328]. 

In Al-doped silicon, Ezhevskii et al. [329] found two EPR 
spectra with monoclinic-I symmetry (NL40 and NL41, see 
Appendix A). A resolved hyperfine structure due to 27 Al was 
observed for some EPR orientations. Although the identifi- 
cation of the centers by using the 57 Fe isotope has not been 
done, the low symmetry of the center and the high value of 
the spin 5 = 5/2 enabled Ezhevskii et al. [329] to identify 
the centers as (Fej*"Fef Al") in two different configurations. 
Hie hyperfine splitting due to two ^Fe iron atoms was later 
reported by Irmscher et al. [215], who also showed that the 
FeFeAl defect produces two groups of EPR lines correspond- 
ing to the ground and the excited states of the same defect 
They reported that storage of the samples at room tempera- 
ture led to the disappearance of the EPR signal, which could 
be restored by shining above-band-gap light on the samples. 



Thus, the EPR studies revealed at least 8 complexes (see 
Appendix A) consisting of two or four iron atoms. To the best 
of our knowledge, there are no experimental data on electrical 
and recombination properties of these centers. 



6 Summarizing discussion 

Research on properties of iron in silicon has almost as long 
a history as the studies of silicon itself. Iron was a com- 
mon unintentional contaminant in the first investigations con- 
ducted in the 1960s/ 1970s. Many experimentalists reported 
the electrical levels of "quenched-in defects" that were ac- 
tually the electrical levels of iron and its complexes. In 
the beginning of the 1980s it was understood that there 
are three major sources of unintentional iron contamina- 
tion: iron contamination during silicon growth (see, for ex- 
ample, [160,330]), surface contamination of wafers during 
chemical cleaning and handling (see, for example, [331]), 
and introduction of iron from the quartz tubes and ambient 
during high-temperature anneals (see, for example, [154, 156, 
332,333]). Unintentional iron contamination levels after heat 
treatment were reported to be as high as 5 x 10 15 cm" 3 [73, 
96,157,160,334]. Later research showed that these contam- 
ination levels can be easily avoided by using clean room 
facilities, high purity chemicals, and state-of-the-art silicon 
wafers. 

A number of studies of iron intentionally introduced in sil- 
icon were done between 1980 and approximately 1990. How- 
ever, the cleanliness of the sample preparation and the quality 
of the as-grown material was not yet sufficient to exclude the 
possibility of cross-contamination with other impurities (pos- 
sibly grown-in to the wafers) which can form complexes with 
iron. In the beginning of 1990s it has become possible to keep 
the level of undesired impurities at the level of 10 n cm~ 3 and 
below, which makes it possible to refine the data on iron in 
silicon accumulated during the past 30 years, and to perform 
new experiments that were not feasible in the past. 

Investigations of iron-related defects provide a perfect ex- 
ample of how different experimental techniques, sensitive to 
structural, recombination, electrical and optical properties of 
defects, can be combined to identify the nature of each de- 
fect. These techniques include studying samples with differ- 
ent types and levels of doping (variation of the Fermi level 
position and charge state of the defect), establishing corre- 
lations between the amplitude of the signals and the anneal- 
ing behavior, correlating ionization energies, illuminating the 
samples (kinetics of dissociation of pairs or switch to another 
metastable state), comparing experimental results with theor- 
etical calculations, and doping silicon with isotopes of iron 
and other impurities under investigation. 

Research over the years showed that the quenching rate of 
the samples is enormously important and often determines the 
type of complexes which are formed in the sample [79,335, 
336]. For example, Wiinstel and Wagner [79] showed that 
after a quench of a p-Si sample in silicon oil, only two iron- 
related peaks were observed (Ey + 0. 1 eV and Ey + 0.43 eV). 
A slower cool (air cool) resulted in the appearance of three 
new peaks at Ey + 0.33 eV, Ey + 0.4 eV, and at Ey + 0.52 eV. 
Finally, after a furnace cool, only a peak corresponding to 
Ey + 0.33 eV was detected. Many iron-related defects require 
special activation or treatment of the sample before they can 
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be detected. For instance, some of the defects appear only in 
irradiated silicon or require illumination to become activated, 
some increase in density during anneals at slightly elevated 
temperatures, whereas the others are unstable and dissociate 
(or grow into larger agglomerates) within hours after diffu- 
sion of iron and quench to room temperature. 

This review was intended as a summary of the past funda- 
mental research of iron in silicon and as a starting point for 
future investigations. We showed that significant progress has 
been made since the 1983 review of one of the authors [19]. 
Appendices A and B provide a comprehensive list of EPR 
resonances and deep levels associated with iron and its com- 
plexes. However, many problems are still unresolved. The 
current knowledge of iron in silicon can be summarized as 
follows. The electrical properties of interstitial iron are well 
understood. The donor level of iron at Ey + 0.38 eV as meas- 
ured by Hall effect and DLTS is well established, and capture 
cross-sections for holes and electrons are available (see Ap- 
pendix B). However, the temperature behavior of the iron 
level in the band gap requires further investigations, particu- 
larly at temperatures below 800 °C. Since the energy position 
of iron at elevated temperatures determines its solubility and 
effective diffusivity in p+ substrates, and thus the segregation 
coefficient in p/p+ gettering, the temperature dependence of 
the position of the iron level in the band gap is not only of 
fundamental interest, but is also highly relevant for industrial 
processes. 

The photoluminescence signal at a photon energy of 
735 meV («1.7ujn) [337-339] and initially identified as 
that of interstitial iron needs verification. Singh et al. [337] 
and Weber et al. [118] estimated the energy position of the 
ground state of the defect responsible for the 735 meV lu- 
minescence band at E v + (0.39 to 0.44) eV, which coincides 
with the position of the energy level of the interstitial iron. 
However, Schlesinger et al. [340] pointed out that there is no 
detectable isotope shift in the previously identified "phonon 
replicas", which they expected for interstitial iron. 

The solubility of iron in intrinsic silicon is well estab- 
lished. However, the data on iron solubility in moderately and 
heavily doped /i-type silicon are limited to a few doping levels 
and temperatures, and need further investigations. The ex- 
perimental data on iron solubility in p-Si and p + -Si are also 
limited. The existing theory of solubility enhancement due to 
the Fermi level effects and pairing with shallow acceptors can 
in principle be used to calculate the solubility for any given 
temperature or doping, provided the temperature dependence 
of the iron level in the band gap is well established. 

The data on the diffusion coefficient of iron are numerous 
and, with a few exceptions, agree with each other. However, 
the existing expression for the diffusivity of iron does not dis- 
tinguish the diffusivity of iron in ionized or neutral charge 
states, and represents the "apparent", or "effective" diffusion 
coefficient There are experimental indications and theoretical 
predictions that the diffusion coefficient of iron does depend 
on the charge state. However, the data where the charge state 
is explicitly specified are obtained in narrow temperature in- 
tervals, are strongly scattered, and often contradict each other, 
as discussed in Sect. 1.3.2. Thus, no conclusion on diffusiv- 
ity of iron in ionized and neutral charge state can be done at 
the present state of knowledge. Furthermore, the dependence 
of the iron diffusivity on the doping level of the substrates has 
never been systematically studied. 



Substitutional iron has never been detected by EPR or 
electrical measurement techniques such as Hall effect or 
DLTS. There are numerous hypothetical assignments of par- 
ticular electrical levels to Fes, but none of them is convinc- 
ingly supported by experimental data. It is well known from 
theoretical considerations [36] that the equilibrium concen- 
tration of Fes should increase with increasing n-type dop- 
ing, if Fes forms acceptor levels (S) in the band gap. It was 
suggested that the experimentally observed enhancement of 
iron solubility in n + -Si is due to the substitutional fraction 
of iron [35]. However, there are not enough studies of n + 
silicon to make any definite conclusions. The only experi- 
mental technique that claims to detect Fes is MSssbauer spec- 
troscopy (see Sect. 1.1), which uses ion implantation that may 
be subsequently responsible for an enhanced substitutional 
fraction of iron. 

Pairing of iron with shallow acceptors, particularly with 
boron, was addressed in many publications and substantial 
progress has been made in both experimental and theoret- 
ical areas of the studies. The process of formation and dis- 
sociation of pairs is well understood qualitatively. However, 
a quantitative description of the pairing kinetics is hindered 
by significant disagreement between various studies in the ac- 
tivation energies for the rate of pairing and dissociation (see 
Sect. 2). The donor energy levels of the pairs in the lower half 
of the band gap are well established, with the exception of 
the level of the orthorhombic configuration of the FeB pair 
which still needs to be confirmed. The acceptor levels of the 
iron-acceptor pairs in the upper half of the band gap are only 
tentatively assigned. 

An important issue that has already been addressed by 
several theoreticians is the type of bonding in the pairs. Al- 
though there is no doubt that the electrostatic interaction 
determines the pairing kinetics, there are indications that co- 
valent or mixed type of bonds are established once a pair is 
formed. However, the exact role of the covalency is uncertain 
since conclusions of different groups contradict each other. 
The scope of the models ranges from purely ionic bonding, 
covalent bonding via the surrounding silicon atoms, to direct 
covaleht bonding between the iron and the acceptor. Appar- 
ently, more detailed theoretical studies using better models 
and larger clusters are required. The requirement of larger 
clusters follows from the ENDOR studies which indicated 
that the iron interacts with up to 100 of the surrounding sili- 
con atoms [24,53,58]. 

There was an intensive discussion in the literature about 
the identification of the photoluminescence signal at 1067 
meV ( ^ 1.16 ujn) as originating from iron and boron [341]. 
Since this spectrum was found exclusively in boron-doped 
silicon, which was intentionally or unintentionally doped 
with iron, it was initially interpreted as the signal of FeB 
pairs. However, a number of experimental results obtained 
later contradict this simple interpretation. First, there is no 
definite correlation of the intensity of the PL signal with 
the concentration of iron or boron [341]. Further, Zee- 
man splitting was nearly isotropic and exhibited no indica- 
tion of an axial nature, characteristic for FeB pairs [341], 
Schlesinger and McGill [342] reported that there was no 
isotopic shift in the no-phonon line or in the phonon repli- 
cas after diffusion of the isotopes ^Fe and ^Fe, A shift 
would be expected if the phonon mode was a local vi- 
bration mode of iron as was suggested by Sauer and We- 
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ber [244]. In accordance, Kluge et aL [343] reported that 
anneals of the quenched boron-doped samples at 60 °C and 
at 90 °C resulted in an increase in the intensity of the FeB- 
related EPR signal, whereas the intensity of the PL signal 
at 1.16p,m decreased by more than an order of magni- 
tude. They suggested that the 1.16 jun PL signal is not due 
to FeB pairs, but rather due to complexes which include 
several iron or boron atoms. This discussion was recently 
resolved by Henry et al. [344], who convincingly showed 
that this PL line is in fact related to gold. Thus, it is pos- 
sible that this line is not related to iron at all and the whole 
entirety of the FeB-related photoluminescence studies is 
just a lesson in the importance of controlling unintentional 
contamination. 

The same ambiguities apply to the identification of the 
PL signal at around 1.1 18 eV, tentatively associated with iron 
and indium [240-242,247,345-348]. Weber et al. [244,347] 
concluded from Zeeman measurements that the PL center in 
In-doped silicon consists of at least two different atoms and 
has (100) symmetry, consistent with the symmetry of Fein 
pair. Schlesinger et al. [248] reported that the PL signal in 
Si:In increases with increasing Fe concentration. In contrast, 
Watkins et al. [243] failed to find any correlation of the PL 
intensity with the iron concentration. 

The viability of the pairing of iron with phosphorus, and 
with shallow donors in general, remains a puzzle. There are 
data that, on one hand, suggest that such pairs do form and 
are stable, and on the other hand, indications that there is no 
pairing whatsoever. This problem will be probably addressed 
more carefully in the near future since the pairing of transi- 
tion metals with shallow donors is an important mechanism 
in the explanation and modeling of phosphorus diffusion get- 
tering of iron [349-353], one of the most efficient and least 
understood gettering techniques. 

The high reactivity of iron results in the formation of 
complexes with other impurities as well. However, only com- 
plexes of iron with gold have been sufficiently studied experi- 
mentally, while the feasibility of pairing of iron with platinum 
has been studied only theoretically. Considering the impor- 
tance of iron and oxygen in silicon, it is surprizing how lit- 
tle research has been done on the interaction between iron 
and oxygen. It seems that both the complexing of iron with 
interstitial oxygen and the influence of iron on the nucle- 
ation and growth of oxygen precipitates are feasible. How- 
ever, there exist only preliminary experiments and hypotheses 
rather than a clear understanding of these processes. The in- 
teraction of iron with carbon, the second major impurity in 
CZ silicon, has never been systematically studied at all (we 
are aware of only one paper, discussed in Sect. 4.4). 

Hydrogen passivation of grown-in defects is used in the 
industry to improve the quality of solar-grade silicon and is of 
significant importance for fundamental research. Complexing 
of metals with hydrogen may result in either the disappear- 
ance of existing levels, or in the formation of new defect 
levels associated with the metal-hydrogen complexes. The 
length of the list of iron-related defect levels in Appendix B 
suggests that a complete study of the interaction of hydrogen 
with each of the iron-related complexes will be a complicated 
and a time-consuming task. Currently, only few experiments 
have been reported. There are indications that hydrogen does 
not passivate interstitial iron, but passivates complexes of iron 
(see Sect. 4.5). 



EPR studies discussed in Sect. 5 demonstrated that iron 
forms complexes consisting of several iron atoms, such as 
FeiFetV, FeiFeiB, Fe^Al, FejFej, and (Fei) 4 . One can ex- 
pect that at least a substantial part of iron agglomerates listed 
in Appendix A form deep levels in the silicon band gap. Stud- 
ies of the electrical activity of these levels have not yet been 
done. Some of these complexes, as well as pairs of iron with 
other impurities may have already been detected in electrical 
measurements, although are not yet identified. They are pri- 
marily deep levels in the upper half and in the middle of the 
band gap. A number of researches reported a level in iron- 
contaminated samples at the middle of the bandgap, around 
£ v + (0.55 ± 0.05) eV, or E c - (0.55 ±0.05) eV. This level 
was originally reported by Collins and Carlson [26] and then 
discussed in [52,79,82,87,88,124,335,354-361]. The con- 
centration of the level was about 1% of the iron concentration 
determined by the diffusion temperature [87,355]. Since the 
energy position of the mid-gap iron-related level reported in 
different references scatter by up to 0.1 eV, it is difficult to 
say whether the same level was observed in all reports. It 
is interesting to note that Wunstel and Wagner [79] detected 
this peak only after a relatively slow air cool of the sam- 
ple, but failed to detect it after both fast cool (oil quench) or 
very slow cool (cool in the furnace). Furthermore, they ob- 
served the level at Ey + 0.52 eV only in boron-doped silicon 
and not in gallium-doped silicon. This led them to conclude 
that boron was involved in the formation of the Ey + 0.52 eV 
level, Thus a possible candidate for this center is one of 
the FeFeB complexes, detected by EPR (see Appendix A). 
Muller et al. [29] speculated that the mid-gap level may be 
due to clusters of several iron atoms, i.e., (Fej) 2 or (Fe04. 
Several other unidentified levels that are reportedly associ- 
ated with iron are located in the upper half of the band gap 
(see [95,161,362-367]). 

Thus, while very much is already known about iron in 
silicon, it will definitely remain an exciting topic for funda- 
mental and applied physical studies for the years to come. 
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Appendix A 

A Table of EPR parameters of iron-related defects in silicon. 
The orientation of the principal axes of the g-tensor is as fol- 
lows; g { is parallel to the (100) axis, g 2 and g$ are parallel 
to (011) and {Oil) axes, respectively. An additional angle, 
6\ determining the angle of rotation of the coordinate system 
around <011) axis is given for the complexes with monoclinic 
symmetry. 
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Defect 


Symmetry 


Spin S 


5i tea) 


Principal g values 
52 (5±) 


53 (5i) 


Reference 


Fe? 


cubic 


1 


2.0702 


2.0702 


2.0702 


[52] 


(a) 






2.0699 


2.0699 


2.0699 


[28,51,56] 


Fef 


cubic 


3/2 


3.524 


3.524 


3.524 


[28,56] 


(b) 














(FejV) + 


trigonal 


3/2 


2.1163 


2.0935 


2.0935 


[29] 


(NL19) 














(d) 














(FejBO 0 


trigonal 


1/2 


2.0676 


4.0904 


4.0904 


[168,239] 


(c) 


(D = 2.7 cm" 1 ) 


3/2 


2.0676 


2.0452 


2.0452 


[164,168] 






3/2 


Z068 


2.0 


2.0 


[28] 


(FeiBs)- 


trigonal 


1 


2.1345 


2.1345 


2.1345 


[167] 


(c,e) 


(£>= 1.38 cm" 1 ) 


1 


Z080 


2.070 


2.070 


[168,216] 






1/2 


4.090 


2.068 


2.068 


[166] 


(FejAls)"** 


trigonal 


1 


4.205 


2.3 


2.3 


[168,215,216] 


(c) 


(Z) = -23cm- 1 ) 


1/2 


8.395 


1.076 


1.076 


[217,218] 


(FejAIs) 0 


trigonal 


1/2 


6.389 


1.138 


1.138 


[213] 


(NL27) 




1/2 


6.366 


1.144 


1.144 


[168,215] 


fee) 




1/2 


6.389 


1.138 


1.138 


[239] 






1/2 


6.370 


1.142 


1.142 


[166] 




fri a final 




2.075 


2.075 


2.075 




fee 1 ! 


(D — 1 46cm~h 












(FeiAls) 0 


orthorhombic-I: 












(c) 


ground state 


1/2 


5.885 


1.236 


1.612 


[213,214,239] 




/XTT 00\ 














fxH tprl ctate 


1/2 


5.877 


1.242 


1.614 


ri66i 

Liooj 




description 














of both lines 


1/2 


1,73 


2.51 


5.36 


[214,239] 




(D = -6.5 cm-', 














E/D = -0.2540) 


3/2 


2.080 


2.059 


2.043 


[214,215] 




trigonal 


1/2 


5.087 


2,530 


2.530 


r 9» 9501 






1 /9 


J.U07 


Z-J-jU 




|PO,Zl't| 


(FejGas) 0 


orthorhombic-I: 












(c) 


ground state 


1/2 


6,19 


0.59 


0.69 


[214,239] 






1/2 


2.02 


4.65 


j-j i 


[214,239] 




Hpc/vintinn 














of both lines: 














(£/D = -0.09, 


3/2 


Z07 


1.97 


105 


[214] 




D« -6 cm" 1 ) 
















1 /I 


£ 38 


1 OR 


1 OR 

l.UO 


r997— 9V11 


(c.f) 














(Lu2) 














(Fejlns) 0 


orthorhombic I: 












(c) 


ground state: 


1/2 


6.26 


0.36 


0.35 


[226,230] 




excited state: 


1/2 


Z070 


3.78 


4.40 


[28,239] 






1/2 


r\ /V7A 

2.070 


3.80 


4.42 


r*itA io£ 01m 
[214,22o,230J 




UCSvl ipUUll 














of both lines 


3/2 


2.09 


2.07 


2.05 


[227] 




(£> = -4.9 cm- 1 , 














£/£> = 0.052) 












(Feilns) + 


orthorhombic 


1/2 


<0.3 


14.1 


<0.3 


[227,232] 


(c) 














tentatively 














FejAus 


trigonal 


1/2 


2.0993 


2.1165 


2.1165 


[278] 






1/2 


2.1001 


2.1178 


2.1178 


[272] 




ICIlutUVwiy 


1 /9 


9 mn 


9 t9fi 

Z. 1 i\> 




f9GAI 
tz^OJ 


(tentatively 


monoclinic-I 












(FeS)") 


0 = 15° 












reps 


tentatively 


1 /9 
If I 


1 019 

1 .7JO 


Z.UI J 


9 OAT 
Z.yOZ 


12V0J 


(tentatively 


monoclinic-I 












(FeS) 0 ) 


0 = 28.5° 












(a) 














FeiSsOH) 


tentatively 


1/2 


1.991 


2.042 


1503 


[296] 


(tentatively 


monoctinic-I 












(FeS) + ) 


0 = 46° 
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Defect 


Symmetry 


Spin 5 




Principal g values 


53 (£J.) 


Reference 


FejSs(IV), 


tentatively 


1/2 


1.939 


1.956 


2.691 


[297] 


possibly 


orthorhombic 












Fei(S-S) 














(c) 














FePd 


trigonal 


1/2 


2.0407 


2.0887 


2.0887 


[57,289] 


FePt 


trigonal 


1/2 


2.0124 


2.1264 


2.1264 


[57,289] 


Fe-oxygen 


trigonal 


1/2 


1.9962 


1.9956 


1.9956 


[309] 


(g) 












Fe-oxygen 


orthorhombic-I 


? 


1.9928 


1.9962 


1.9926 


[309] 


tt> 












FeOV 


monoclinic-l 


1/2 


4.20 


2.15 


4.10 


[57] 




0 is not reported 








(FeiFeiV)" 


trigonal 


1/2 


2.059 


6.235 


6.235 


[291 


(d) 




3/2 


2.059 


2.078 


2.078 


[29] 


(NL20) 










(FtSjFejV) + 


orthorhombic II 


5/2 


2.51 


1.47 


0.57 


f291 


(d) 










(NL25) 














FeiFfeiVa 


monoclinic-I 


1/2 


7.38 


1.961 


4.90 


[29] 




n oo Qo 

u — JJ.O 


< /9 
5/2 


9 t nn 




z.uuo 


[29] 


(NL2I) 










(FeiFei) + 


monoclinic-I 












(n) 


ground state: 


1/2 


1.15 


9.44 


2.06 


[29,223,224,326] 


(NL24) 


CXCllGCL SLoiC. 

6 = 36.4° 


1 /9 


i on 


o.DU 


j.U/ 


[223,224] 




description 


5/2 


«2.05 


^2.05 


5^2.05 


[223] 




of both lines: 












0.20) 














UlgUIloI 






9 /Vi2 
Z.UOo 


9 AAQ 
Z.UOO 


L29] 


/KIT OOX 












FeFeB 


monoclinic-I 












(c) 


ground state (?) 


1/2 


1.49 


2.91 


8.97 


[223,224] 










o one 




[328] 






1 /9 








TOO 7 0041 




(7 — JO. / 














hnth linpc* 

UvLU LUIGo. 








e« 9 07 


















FeFeB 


orthorhombic-I: 












(c) 


ground state (?) 


1/2 


7.55 


4.18 


1.81 


[223,224] 


(NL32) 






7 009 


A 1(14 
*r. lOf 


1 ft 1 1 
1,01 1 






CALllvU aUli-C \l) 


1 /9 


1 7^ 


1 O/L 








H<*c/"*t"i ration nf 


S/9 


as 9 OS 






f90^ "5941 




\\f\tYt liriAc* 

UUUl JltlCo. 














\Jz,f LJ ~ U.Uo J 












FeFeAl 


monoclinic-I 


1/2 


1.470 


2.836 


8.907 


[329] 


(c) 


ground state (?) 


1/2 


1.48 


2.82 


8.89 


[215] 


(NL40) 




1/2 


1.404 


2.661 


9.027 


[329] 




excited state (?) 


1/2 


5.48 


3.17 


2.86 


12151 




0 = 36.6° 












description of 














both lines: 


5/2 


2.064 


2.064 


2.064 


[329] 




(E/Dk 0.1447, 


5/2 


2.081 


2.023 


2.083 


[215] 




£> = 3cm~ 1 ) 













(a) observed predominantly in n-Si. 

(b) observed predominantly in />-Si. 

(c) appears after one to thirty days of storage at room temperature, necessary for association of the pairs. 

(d) observed in irradiated silicon. 

(e) requires activation by above-band-gap light 

(f) activated by illumination with light with hv > 0.5 eV. 

(g) was detected in CZ-grown silicon by EDMR. 

(h) unstable at slightly elevated temperatures (T> 60 °C). 

Appendix B 

Table of electrical levels associated with iron and its com- 
plexes. Energy level positions were determined using deep 



level transient spectroscopy (DLTS), Hall effect (Hall), pho- 
tocapacitance technique (PhCap), electron paramagnetic res- 
onance (BPR) and photo-EPR, Fourier transform infrared 
spectroscopy (FTIR), and thermally stimulated capacitance 
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(TSCAP). Uncertainty in the energy position reflects only 
scatter in experimental data reported by different groups, but 
does not include possible uncertainty of their data (which 
is usually not reported). Capture cross-sections were deter- 



mined from temperature dependence of emission rate (Arrhe- 
nius plot) with exception of several studies that used direct 
isothermal methods (see Chapt 1 1 of [369] for details of the 
experimental techniques). These studies are indicated by (*). 



Defect 


Energy position /eV 
measurement methods 


Hole capture 
cross-section 


Electron capture 
cross-section 


Fe? /+ 


£V + (0.38±O.Ol)eV 
[26,75,78,80-82,88-90] 
DLTS, Hall, PhCap, EPR 


3.9x 10- |6 exp(-0.045eV/Jfc B r) cm 2 
[66,80,90-93] 
(*) 


cr„s»4xl0" l4 cm 2 
(r = 300K) [106,107] 


(FeB)°/+ 
trigonal 


£ v +0.1 eV, 
[80,91,92,124] 
DLTS 


Op*2xl0- 14 cm 2 [80] 
or p «1.5xlO~ 13 cm 2 
(r = 55K)[92] 
Op > *fx iu cm 


<r n «4xl0- ,3 cm 2 [107] 


(FeB)-/°. 
trigonal 


£ c - (0.26 ± 0.03) eV 
[90,91,167,172-174] 
DLTS, Hall, photoEPR, 
FTIR 


tr p as3xl0- 14 cm 2 
(7* = 300 K) [175, 176] 
Op«3x 10" 15 cm 2 
(r = 90K)[90] (*) 


an^2.5xl0- ,3 cm 2 
(T = 300 K) [175,176] 
cr n 10~ 15 cm 2 
(r = 90K)[90] (*) 


(FeB)°/+ 

orthorhombic 

(tentatively) 


£ v + 0.07eV 
[86,186] 
Laplace Transform DLTS 


? 


? 


(FeB)-/°, 

orthorhombic 

(tentatively) 


£ c -0.43eV 
[172,188] 
TSCAP 


? 


ff n ^3xl0- |5 cm 2 , [188] 


(FeAl)°/ + , 
trigonal 


E v + (0.20 ± 0.01) eV 
[78,80,81, 144, 154, 170,223] 
DLTS 


cr p ss2xl0- 14 cm 2 [220] 
Op«4xl0- 14 cm 2 [80] 


? 


(FeAJ)°/ + , 
orthorhombic 


£ v +0.13eV 
[78,80,144,220,223] 
DLTS 


Op « 6 x 10~ 14 cm 2 [220] 


? 


(FeGa)°/+ 
trigonal 


E W + 0.24 eV [80, 170,225] 
DLTS 


a p ^2xl0- 14 cm 2 [80] 
a p 3 x 10" 15 cm 2 (T = 120 K) [170] 


? 


(FeGa)°/+ 
orthorhombic 


Ev +0.14 eV [225] 
DLTS 


? 


? 


(Feln) 0 /+ 
trigonal 


K v +0.27eV [2251 
DLTS 


? 


? 


(Fein)-/ 0 , 
trigonal 


£ c -0.32eV [233] 
FTIR 


? 


? 


(Felr0<>/+, 
orthorhombic 
(Fein)"/ 0 , 
orthorhombic 


£ v +0.I5eV [225] 

DLTS 
Eq — 0.39 ev [233] 

FTIR 


? 
7 


? 

*? 


(FeAu) u/ ^ 


Ey + (0.44 ± 0.01 ) ev [157,2oO,2o2J 
DLTS, Hall 


a p 2 x 10 cm [282] 




(FeAu) _/0 


£ c -(0.34db0.01)eV [280,283,284] 
DLTS 


? 


o n » 3 x 10~ 15 cm 2 [284] 
cr„«2xl0- |5 cm 2 [280] (*) 


FePd 


E c - 0.32 eV [291,292] 
DLTS 


? 


? 


FePd 

(tentatively) 


£ v +0.52eV [292] 
DLTS 


? 


? 


(FeZn)°/ + 


Ey + 0.17 eV [295] 


op = 5 x 10- 15 cm 2 (T = 90 K), [295] 


? 


(FeZn)-/° 


£ c -0.47eV [295] 


? 


a n = 1.5x 10- |6 cm 2 
(T = 250K),[295] 


FeOV 


E c - 0.36 eV [304] 
DLTS 


? 


ff n «4xl0- l5 cm 2 [304] 


Fe-0 
(tentatively) 


£ v +0.33eV [79] 
DLTS 


? 


' ? 
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